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Avant-propos  . 


Monsieur  le  Professeur  Quentin  et  les  organisateurs  de  ce  colloque  m'ont 
accorde  le  grand  honneur  de  presenter  le  premier  expose  scientifique  de  ce 
Colloque  sur  la  Diffusion  des  Ondes  Ultrasonores ,  ainsi  que  les  exposes 
introductifs  de  chaque  partie  du  colloque.  j'apprecie  amplement  cet  honneur  qui 
me  permet  de  presenter,  dans  le  cadre  de  la  venerable  Ecole  Normale  Superieure 
de  l'Universite  de  Paris,  a  une  a  semblee  de  chercheurs  renommes  et  d'etudiants 
qui  le  seront  un  jour  egalement,  les  fruits  d'un  travail  de  vingt  ans  sur  les 
problemes  d'acoustique  que  j'ai  pu  effectuer  avec  mes  collaborateurs  et  etudiants 
au  Departement  de  Physique  de  l'Universite  Catholique  d'Amerique  a  Washington, 
ainsi  qu'avec  de  nombreux  chercheurs  dans  les  laboratoires  gouvernementaux  des 
environs  de  Washington  avec  qui  j'ai  eu  la  chance  d' avoir  pu  col laborer. 

Dans  la  plupart  de  mes  exposes,  je  me  suis  donne  la  peine  de  non  seulement 
presenter  un  resume  de  ce  qu'on  a  fait  anterieurement,  mais,  puisque  nous  sommes 
toujours  en  plein  travail,  d'y  inclure  aussi  les  derniers  r£sultats  de  nos 
recherches  qui  ne  sont  pas  encore  publies,  afin  de  mettre  au  courant  les 
participants  du  colloque  de  nos  interets  actuels.  D'une  grande  importance,  aussi, 
est  l' inclusion  dans  mes  exposes  de  quelques  resultats  de  nombreuses 
investigations  experimentales  faites,  ici,  dans  les  laboratoires  Franqais, 
investigations  qui  ont  ete  inspirees  par  nos  theories  anterieures,  et  qui 
heureusement  les  ont  bien  verifiees,  telles  les  experienres  faites  par  le 
Professeur  Quentin,  le  Docteur  de  Billy  et  leur  equipe  a  l'Universite  Paris  VII 
sur  la  penetration  d'objets  elastiques  par  l'ultrason  de  haute  frequence  ;  et 
les  investigations  sur  les  resonances  et  les  ondes  de  surface  sur  des  objets 
elastiques  faites  par  le  Professeur  Ripoche,  le  Docteur  Maze  et  leurs  collegues 
a  l’Universite  du  Havre,  soutenus  theoriquement  par  les  Docteurs  Derem  et 
Rousselot  du  Laboratoire  Central  des  Telecommunications  a  Velizy-Villacoublay  ; 
plus  recenment  aussi  les  experiences  des  Docteurs  Gazanhes  et  Sessarego  au  CNRS 
de  Marseille,  et  autres  recherches. 

Nous  sommes  profondement  satis faits  que  nos  suggestions  theoriques  aient 
cause,  si  je  peux  dire  ainsi,  autant  d'echos  en  France.  J'espere  tout 
specialement  qu'a  partir  de  nos  etudes  sur  les  resonances,  et  les  experiences 
associees,  de  nouvelles  methodes  et  applications  vont  pouvoir  etre 
developpees  dans  le  domaine  acoustique. 


Herbert  Uberall 

Paris,  le  21  Novembre  1984 


Pre  face 


Professor  Quentin  and  the  organizers  of  this  colloquium  have  honored 
me  by  their  request  that  I  present  the  first  scientific  talk  in  this 
Colloquium  on  the  Scattering  of  Ultrasonic  Waves,  as  well  as  the  opening 
talks  of  each  session  of  the  colloquium.  I  greatly  appreciate  this  honor 
which  permits  me  to  present,  in  the  surroundings  of  the  venerable  Ecole 
Norma le  Sup^rieure  of  the  University  of  Paris,  to  an  assembly  of  renowed 
scientists  and  of  students  aspiring  to  such  a  status  also,  the  fruits  of 
twenty  years'  work  on  acoustic  problems  which  I  have  been  able  to  carry 
out  with  my  collaborators  and  students  at  the  Department  of  Physics  of 
the  Catholic  University  of  America  in  Washington  DC,  as  well  as  with 
numerous  scientists  in  Government  laboratories  of  the  Washington  area, 
with  whom  I  had  the  good  fortune  to  be  able  to  collaborate.  In  most  of  my 
lectures,  I  took  pains  to  present  not  only  a  resume  of  what  was  done  by  us 
previously,  but,  since  we  always  keep  on  working,  to  include  also  the 
latest  results  of  our  studies  which  have  not  yet  been  published,  in  order 
to  keep  the  participants  of  the  colloquium  informed  of  our  present  interests. 
Of  great  importance,  additionally,  is  the  inclusion  in  my  talks  of  some  of 
the  results  of  numerous  experimental  studies  that  were  carried  out  here,  in 
French  laboratories,  studies  which  were  induced  by  our  earlier  theories,  and 
which  they  fortunately  have  well  verified,  such  as  the  experiments  of 
Professor  Quentin,  Dr.  de  Billy  and  their  team  at  the  University  of  Paris  VII 
on  the  penetration  of  elastic  objects  by  high-frequency  ultrasound;  and  the 
experiments  on  resonances  and  surface  waves  on  elastic  objects  carried  out  by 
Professor  Ripoche,  Dr.  Maze  and  their  colleagues  at  the  University  of 
Le  Havre,  supported  theoretically  by  Drs.  Derem  and  Rousselot  of  the 
Laboratoire  Central  des  Telecommunications  of  Veiizy-Villacoublay^  and 
recently  also  the  experiments  of  Drs.  Gazanhes  and  Sessarego  at  the  CNRS  at 
Marseilles,  and  those  of  other  scientists.  We  are  profoundly  satisfied 


that  our  theoretical  suggestions  have  caused,  if  I  may  say  so,  such  an 
echo  in  France.  I  specially  hope  that  out  of  our  studies  on  resonances, 
and  of  the  related  experiments,  new  methods  and  applications  will  be 
developed  in  the  field  of  Acoustics. 
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RESUME  :  On  obtient  l'expression  du  faisceau  acoustique  borne  le  plus  general, 

c'est-S-dire  celle  d'un  faisceau  fini  dans  deux  dimensions,  emergeai  t  d'une 
ouverture  collimatrice  de  forme  arbitraire,  et  possedant  un  profil 
arbitraire  d ' amplitude.  On  etudie  la  divergence  geometrique  d'un  tel  faisceau 
pour  deux  cas  de  profil  d'amplitude  :  un  profil  gaussien  et  un  profil  plan- 
rectangulaire.  Dans  chaque  cas,  on  obtient  la  divergence  geanetrique  oomme  fonctio 
du  rayon  correspondant  ala  dimension  de  1 'ouver-ture  collimatrice  et  de  trois 
parametres  differents  :  longueur  d'onde,  distance  axiale  du  plan  de 
l'ouverture,  et  distance  radiale  de  l'axe  du  faisceau.  Comme  application 

de  cette  theorie  gendrale,  on  obtient  l'amplitude  de  diffusion  d'un 
faisceau  borne  par  une  sphere  rigide  d'une  fa^on  exacte  pour  les  deux 
profils  mentionnes  du  faisceau;  et  on  fait  une  comparaison  entre  les 
resultats  concernant  la  section  efficace,  et  ceux  obtenus  par  la  methode  de 
Kirchhoff.  Mention  est  faite  d'autres  representations  d'un  faisceau  borne 
rencontrees  dans  la  litterature. 

SUMMARY  :  We  derive  the  expression  for  the  most  general  bounded  acoustic  beam, 

i.e.,  for  a  beam  finite  in  two  dimensions,  emerging  from  an  arbitrarily 
shaped  collimating  aperture  at  an  arbitrary  angle  and  with  an  arbitrary 
amplitude  profile.  We  investigate  the  spreading  of  such  a  beam  for  two 
amplitude  profiles  :  a  Gaussian  profile  and  a  flat-topped  square  profile. 

In  both  cases  the  spreading  is  determined  as  a  function  of  the  ratios 
between  the  collimating  aperture  sizes  and  three  parameters  : 
wavelength,  axial  distance  from  the  aperture  plane,  and  radial  distance 
from  the  beam  axis.  As  a  specific  application  of  this  general  theory,  we 
evaluate  the  scattering  of  a  finite  beam  from  a  rigid  sphere  in  an  exact 
fashion  for  the  two  mentioned  amplitude  profiles  of  the  finite  beam,  and 
compare  results  for  the  bistatic  cross  sections  to  those  obtained  using 
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the  Kirchhoff  approximation.  We  also  mention  other  representations  for 
finite  beams  encountered  in  the  literature. 

1.  INTRODUCTION  : 

Des  faisceaux  bornes  ont  ete  etudies  theoriquement  et  experimentalement. 
Schoch  [1,2]  a  etudie  des  effets  de  deplacement  apparents  d'un  faisceau 
borne  .reflechi  par  un  solide  ;  d'autres  etudes  sur  ce  sujet  ont  ete 
faites  par  Mayer  et  al  [3,4]  et  Breazeale  et  al  [5,6].  Les  experiences  de 
Neubauer  [7],  et  les  theories  de  Tamir  et  Bertoni  [8,9]  ont  clarifie  le 
role  joue  par  l'onde  de  Rayleigh  dans  ces  phenomenes.  Plus  recemment, 
une  autre  theorie  a  ete  faite  par  Claeys  et  Leroy  [10]  concernant 
1 'interaction,  d'un  faisceau  borne  avec  la  surface  plane  d'un  solide  ; 

la  reflexion  d'un  faisceau  bom&  par  une  surface  rugueuse  (periodique) 
a  ete  etudiee  egalement  [11]. 

Dans  tous  ces  travaux,  on  a  considere  des  faisceaux  acoustiques 
bornes  dans  une  dimension  seulement,  ce  qui  etait  suffisant  pour  la 
description  theorique  des  experiences  mentionnees.  Pourtant,  1' extension 
de  la  theorie  au  cas  d'un  faisceau  borne  dans  deux  dimensions  ne  parait 
pas  tres  difficile  ;  les  descriptions  de  tels  faisceaux  ont  ete  donnes 
par  Gaunaurd  et  (iberall  [12]  et,  independamment ,  par  Alais  et  Hennion  [13]. 
Toutes  les  theories  representent  les  faisceaux  bornes  par  une  superposition 
continue  d'ondes  planes,  dans  le  sens  de  Fourier,  se  propageant  dans 
des  directions  diverses.  La  seule  exception  est  faite  par  Claeys  et 
Leroy  [10]  qui  representent  le  faisceau  borne  (a  une  dimension)  par  une  superpos 
tion  discrete  d 'ondes  planes’  inhomogenes  [14,15],  se  propageant  dans  la  meme 
direction.  Cette  representation  a  l'avantage  d'une  description  analytique 
du  faisceau  reflechi,  d'une  faqon  analogue  a  la  methode  de  Prony  [16], 
tandis  que  la  methode  de  Fourier  nous  fournit  des  expressions  integrales. 

Dans  ce  qui  suit,  on  developpe  la  theorie  d'un  faisceau  acoustique 
borne  dans  deux  dimensions  basee  sur  la  methode  de  Fourier,  et  on  obtient  la 
divergence  georaetrioue  d'un  f  ais  ceau  sort  ant  d'une  ouverture  collimatrice .  Comme 
application  de  notre  representation  du  faisceau  borne,  on  considere  la 
diffusion  d'un  tel  faisceau  par  une  sphere  ;  les  resultats  exacts  obtenus 
ainsi  pour  la  section  efficace  de  diffusion  peuvent  etre  compares  avec 
ceux  de  la  diffusion  d'une  onde  plane  par  la  sphere,  et  ceux  de  la 
diffusion  du  faisceau  analysee  par  la  methode  de  Kirchhoff. 

-  2  - 


-■  •  ■  i  --*■  *  '  ---• 


2.  FAISCEAUX  BORNES  DANS  DEUX  DIMENSIONS. 


Un  faisceau  infini. 


poo &)  =  e 


(2.1) 


incident  de  -  »  sur  un  plan  z  =  -  L,  a  pour  expression 


Poo  -  r0  rf )  e  oVe 

2  =  -L 


(2.2) 


dans  une  ouverture  de  ce  plan,  ou 


k„  =  k,»). 


,  z),  -n.  =  f?/-L). 


(2.3a) 


(2.3b) 


On  a  introduit  ici  le  profil  PQ(p)  de  1' amplitude  du  faisceau  a  travers 
I' ouverture. 

Le  faisceau  le  plus  general  peut  maintenant  etre  represente  par 
1' expression  de  Fourier, 


f Y(Z)  e*kr  d1-*, 


(2.4) 


avec  la  relation  de  dispersion 


kx  =  ±  (kz-^<‘-)i/i,  z‘  =  z  +  L  2  0  (>C^  C.5a) 


li  =  i  •i(52-k1)V\  2'-  z  +  L  2o  (^l>kl),(2.5b) 


les  signes  etant  choisis  pour  assurer  la  convergence  quand  |z'|  -*■  00 


L' egression  (2.4)  est  une  solution  de  l'equation  de  Helmholtz,  se 
reduisant  correctement  21  l'onde  dans  l'ouverture,  Eq.(2.2),  si  ¥  (x)  est 
la  transformee  de  Fourier  de  cette  onde. 


-I 


Y  &)-  P0  ( ?  )£ 


d2<o. 


(2.6) 


Les  equations  (2.4)  -  (2.6)  precisent  entierement  le  faisceau  borne  dans 
deux  dimensions,  si  on  considere  l'argument  de  Claeys  et  Leroy  [10]  que 
pour  des  faisceaux  larges,  3Pq/3z  peut  etre  presume  comme  proportionnelle 
£  PQ.  Dans  tous  les  travaux  anterieurs,  les  expressions  de  Fourier  comme 
(2.4),  ou  d'une  superposition  analogue  d'ondes  inhomogenes  [10]  etaient 
considerees  dans  une  dimension  seulement,  et  representaient  done  un 
faisceau  borne  dans  une  dimension,  e'est-a-dire  un  faisceau  en  forme  d'un 
rub an  plan. 

Pour  un  faisceau  a  symetrie  de'  revolution,  P  =  P  (p) ,  se  propageant 

o  o 

parallelement a  l’axe  z  (3Cq  *  0),  l’equation  C2.4)_se  simplifie  en 

co  o° 

p(r)  =  (2.7) 


on  peut  alors  considerer  les  cas  speciaux  de 


T0  (q)  =  H  (  -  ?),  (2.8a) 

OO 

‘R0e'LkL  jf  7t  (k!?3)  X  die  (2.8b) 


d'un  faisceau  de  profil  rectangulaire  de  rayon  Rq,  H  etant  la  fonction  de 
Heaviside,  ou 


Zty- 

-  .  ,  -i'kL  f  -~K'<rVZ  t  ,  N  w 

p(-r)  =  <rze  Je  J0(arp)^  -xdK 


(2.9a) 


(2.9b) 


d'un  faisceau  de  profil  Caussien. 
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3.  DIVERGENCE  GEOMETRIOUE  DU  FAISCEAU. 


Se  propageant  en  dehors  de  l'ouverture,  le  faisceau  diverge  I  cause 
de  la  diffraction.  Les  integrales  des  equations  (2.8b)  ou  (2.9b)  peuvent 
alors  etre  evaluees  [12]  par  la  methode  de  col,  ou  par  integration 
partielle . 

La  figure  3.1  presente  en  variables  sans  dimension, la  divergence  geaaetriqu 
d'un  faisceau  rectangulaire  a  l'ouverture  (v  =  0),  a  la  frequence  s  =  kRQ  =  8.  La  figur 
3.2  dome  la  divergence  peametrique  d'un  faisceau  gaussiena  la  freauence  ko  =  &. 

4.  APPLICATION  :  DIFFUSION  D'UN  FAISCEAU  PAR  UNE  SPHERE  0  2} * 

Le  champ  total  d'une  onde  plane  en  presence  d'une  sphere  de  rayon  a 
est  la  somme  de  l'onde  incidente  p.  et  de  l'onde  diffusee  p„„, 

U1C  rsC’ 


P^,k)-  I n  EE  (4..« 

r  '  t--o  L 


1' amplitude  diffusee  C^(ka)  etant,  dans  le  cas  d'une  sphere  rigide  : 


dg[k&)  =  -  jjz  (ka)/ (k<x). 


(4.1b 


La  section  efficace  bistatique  (ou  distribution  angulaire)  de  diffusion 
est  donnee  par 


acr 

-  =  Aiyn  -r 

dSl  f  -+00 


qui  devient  en  champ  lointain 


P.Vt 


(4.2) 


-  (£)'l£<i  (2M)Zfr0j6) 


Un  faisceau  borne  a  pour  champ  diffuse 


±1 

Tt  L 


1--0 


(4.4) 


ce  qui  donne  en  champ  lointain  dans  le  cas  d'un  faisceau  circulaire  : 


U  dff 

4Xx  d0 


(2^1)C(U)L2  (c«e) 

k“  to  *  e  e 


(4 ,5a) 


oQ 


Cos&k)j<d-K, 


(4.5b) 


"X  =  k  sin  0^.  Pour  un  faisceau  de  profil  rectangulaire  : 

I,  =  1  -  f?e(f)eikLU'rVf, 
f  =  {  i  ^K./L)M'i't  , 


(4 .6a) 


(4.6b) 


en  utilisant  la  methode  du  col. 

La  figure  4.1  presente,  pour  une  sphere  rigide,  d0/d9  dans  le  cas  du 
faisceau  rectangulaire,  pour  L  *  1  et  ka  *  3  (1  gauche)  ou  ka  =  5  (5  droite) , 
pour  le  rayon  du  faisceau  r  =  a  ou  r  =  0.866a  (et  egalement  r  *  400  a  ,  en  haut) . 
Pour  comparaison,  la  figure  4.2  montre  les  memes  resultats  obtenus  par  la 
methode  de  Kirchoff  [12,17],  sauf  pour  le  cas  r  >  a  qui  coincide  avec  le 
cas  r  ■  a. 
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5.  CONCLUSION 


Les  calcuiiS  precedents  montrent  que  : 

(a)  un  faisceau  borne  dans  deux  dimensions  peut  etre  represente  par 
une  integrale  bidimensionnelle  de  Fourier  qui 

(b)  permet  un  calcul  de  la  divergence  geanetriaue  dti  faisceau  anergeant  d'  une 
ouverture. 

* 

(c)  rend  possible  1 ' evaluation  de  la  diffusion  du  faisceau  borne  par  une  cible. 

(d)  admet  l'analyse  de  1' interaction  du  faisceau  borne  en  deux  dimensions 
avec  des  structures  telles  que  le  dioptre  ou  la  plaque,  plan  ou  courbe,  ainsi 
que  le  calcul,  par  exemple,  d'un  faisceau  dans  un  milieu  inhomogene.  Une 
partie  des  resultats  rapportes  est  basee  sur  des  travaux  assistes 

f inancierement  par  l'ONR  et  le  NSWC. 
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Fig.  3.1.  Etendue  diffractive  d'un  faisceau  de  profil  rectangulaire , 

avec  s  =  kR  =8. 
o 

Fig.  3.2.  Etendue  diffractive  d'un  faisceau  de  profil  Gaussien 
(x  =  p/a,  y  *  z'/a),  avec  w  =  ka  =  8. 

Fig.  4.1.  Distributions  angulaires  de  la  diffusion  d'un  faisceau  rond, 
avec  profil  rectangulaire,  d'une  sphere  rigide. 

Fig.  4.2  Resultats  comparables  a  ceux  de  la  figure  3,  obtenus  par 
1' approximation  de  Kirchhoff. 
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:  In  this  report,  we  describe  a  complete  resonance  theory  for  the  acoustic 

transmission  and  reflection  coefficients  of  an  elastic  plate  imbedded  in  a 
fluid  medium,  including  effects  of  plate  viscosity.  The  purpose  of  this 
formulation  is  to  provide  a  direct  means  for  determining  the  material 
parameters  of  the  plate  from  the  measured  acoustic  resonances  of  the  Rayleigh 
and  LamL  waves  in  the  plate  (i.e.,  their  positions  in  frequency  or  angle, 
their  widths  and  their  heigths)  which  are  given  in  our  formalism  by  explicit 
analytic  expressions  that  depend  on  the  material  parameters.  Viscosity  is 
seen  to  manifest  itself  in  a  decrease  of  the  resonance  heights  (especially 
for  the  narrow  shear-type  resonances)  and  in  a  broadening  and  frequency 
dependence  of  their  widths,  which  may  be  used  to  determine  the  frequency- 
dependent  loss  factor  of  the  plate.  This  approach  then  solves  the  inverse 
scattering  problem  for  the  case  of  a  plate.  We  also  consider  the  special  case 
of  a  fluid  layer  imbedded  in  another  fluid.  In  addition,  a  layered  ocean 
floor  is  modeled  by  a  sediment  layer  on  top  of  a  denser  substratum,  and 
overlaid  by  the  water  column.  It  is  shown  for  this  case  again  that  resonances 
in  the  acoustic  reflection  coefficient  are  very  prominent  and  can  be  used 
to  determine  the  properties  of  both  sediment  layer  and  substratum. 

:  Dans  cet  article,  nous  presentons  une  theorie  complete  des  resonances 

dans  les  coefficients  de  transmission  et  de  reflexion  acoustique  par  une 
plaque  elastique  immergee  dans  un  milieu  fluide,  en  tenant  compte  des  ef fe ts  dus  a 
la  viscosite  de  la  plaque.  Le  but  de  ce  formalisme  est  de  fournir  un  moyen 
direct  de  determination  des  parametres  du  materiau  constituant  la 
plaque, par  la  mesure  des  resonances  des  ondes  de  Rayleigh  et  de  Lamb 

(c'est-U-dire,  par  leur  position  en  frequence  ou  en  angle,  leurs 
largeurs  et  leurs  hauteurs).  Ces  quantites  sont  donnees  dans  notre  formalisme 


par  des  expressions  analytiques  explicates  qui  contiennent  les  parametres 
du  materiau.  On  trouve  que  la  viscosite  se  manifeste  par  une  decroissance  des 
hauteurs  des  resonances  (part iculierement  pour  les  resonances  etroites  du 
type  transversal),  et  par  un  elargissement  et  une  dependance  en  frequence  de 
leurs  largeurs.  Cette  methode  Dermet  de  rosoudre  alors  le  probleme  de  diffusion 
inverse  pour  le  cas  d'une  plaque.  Nous  considerons  egalement  le  cas  special 
d'une  couche  fluide  immergee  dans  un  autre  fluide.  En  outre,  nous  representons 
les  couches  du  fond  de  1* ocean  par  une  couche  de  sediment  posee  sur  un  fond 
plus  dense,  et  couverte  par  I'eau  de  1' ocean.  Nous  montrons  que,  pour  ce  cas 
egalement,  les  resonances  dans  le  coefficient  de  reflexion  acoustique 
apparaissent  de  faqon  tres  visible,  et  peuvent  etre  utilisees  pour  la 
determination  des  proprietes  de  la  couche  de  sediment  et  du  fond. 

1.  INTRODUCTION. 

A  comprehensive  theoretical  and  experimental  investigation  by  Schoch  [1] 
has  shown  that  the  process  of  reflection  and  transmission  of  an  acoustic  wave, 
incident  through  a  fluid  on  an  elastic  plate  imbedded  in  that  fluid,  is 
capable  of  exciting  elastic  Rayleigh  and  Lamb-type  waves  [2]  in  the  plate. 

Such  an  excitation  takes  place  when  the  trace  of  the  phase  velocity  of  the 
incident  wave,  as  it  moves  along  the  plate,  coincides  [3]  with  the  phase 
velocity  of  the  plate  waves  ;  it  manifests  itself  in  the  appearance  of 
resonancestM  in  the  acoustic  reflection  (R)  and  transmission  coefficient  (T) 
of  the  plate,  either  as  the  frequency  f  is  varied  at  a  fixed  angle  of 
incidence  0  from  the  normal,  or  as  9  is  varied  at  a  fixed  frequency. 

Mathematically,  such  resonances  in  R  or  T  come  about  by  a  vanishing  of 
the  common  denominator  in  the  fractional  analytic  expressions  for  these  two 
quantities.  Zeros  of  this  denominator,  i.e.,  poles  in  R  and  T,  do  not  occur 
for  real  values  of  both  f  and  8  ;  they  occur,  e.g.,  at  complex  values  of 
f  if  8  is  kept  real  (physical),  or  vice  versa.  Since  experimentally  both  f 
and  8  are  real,  measurements  with  variation  of  f  will  carry  us  successively 
past  the  poles  which  are  located  somewhat  off  the  real  f-axis,  thus 
leading  to  finite  heights  and  widths  of  the  observed  resonance  peaks. 

The  resonance  features  in  R  and  T  can  be  studied  by  isolating  their 
pole  contributions,  that  is,  by  replacing  the  denominator  by  its  first  term 
in  a  Taylor  -  series  expansion  about  each  of  the  poles.  This  was  done  by 

,  ,  ,  H  , 

Fiorito,  Madigosky  and  Uberall  for  the  cases  of  a  fluid  layer  in  another 
fluid,  the  layer  being  either  lossless  [A]  or  lossy  [5],  as  well  as  for  an 
elastic  plate  imbedded  in  a  fluid,  with  the  plate  being  lossless  [6]  or 
lossy  [7].  It  is  possible,  however,  to  go  beyond  this  approach,  and  to 
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represent  R  and  T  by  an  exact  series  representation  as  a  sum  of  resonance 
terms  in  the  sense  of  Mittag-Lef f ler  [8]  ;  this  is  also  akin  to  the  theory 
of  resonances  in  nuclear  reactions,  known  as  the  Breit-Wigner  theory  [9]. 

Such  an  exact  representation  is  discussed  here,  as  given  recently  by  Fiorito, 

,  It 

Madigosky  and  Uberall  [10]  for  the  case  of  a  fluid  layer.  The  merit  of  such 
an  exact  resonance  decomposition  rests  in  the  fact  that  it  admits  a 
description  of  the  resonances  even  for  small  impedance  contrasts  between  the 
layer  and  the  ambient  flu^d. 

The  mentioned  investigations  [4-7,10]  hive  provided  analytical  expressions 
for  the  resonance  positions  (in  f  or  0),  their  heights  and  their  widths. 

These  expressions  depend  on  the  materials  parameters  of  the  plate,  and  the 
measurement  of  these  quantities  thus  furnishes  a  direct  means  for  the 
determination  of  the  material  parameters  by  remote  acoustic  sensing.  In  a  way, 
such  a  procedure  then  corresponds  to  a  solution  of  the  "inverse  scattering 
problem",  where  the  properties  of  the  scatterer  are  determined  from  the  form 
of  the  returned  echo. 

An  experimental  study  of  plate  resonances  has  recently  been  carried  out 
by  Izbicki,  Maze,  and  Ripoche  [11—13],  showing  in  exact  detail  the 
backscattering  resonances  in  the  reflection  coefficient  for  aluminum  plates 
immersed  in  water  as  plotted  vs.  f,  and  featuring  their  positions  in  f,  their 
heights  and  widths. 

Fluid  layers  have  also  been  used  to  model  sediment  layers  of  the  ocean 
bottom,  lying  on  top  of  a  substratum  which  was  also  treated  as  a  fluid  [14], 
Calculated  reflection  coefficients  for  a  layered  ocean  bottom  show  very 
prominent  resonance  features  [15].  These  are  explained  in  the  mentioned  model, 
which  proceeds  with  a  detailed  analysis  to  show  what  types  of  measurements 
have  to  be  carried  out  on  these  resonance  features  in  order  to  completely 
determine  the  geometrical  and  material  properties  of  layer  and  substratum. 
Finally,  a  way  of  analyzing  the  layer  resonances  via  the  use  of  long  sound 
pulses,  which  are  employed  in  order  to  cause  a  ringing  of  the  layer  resonances 
one  at  a  time,  has  also  been  pointed  out  [16]. 

2.  THEORY  OF  LAYER  RESONANCES. 

The  exact  expression  for  the  transmission  and  reflection  coefficient  of 
an  elastic  layer  as  calculated  by  Schoch  [1]  is  first  rewritten  in  a 
convenient  form  : 

(  i  1 

=  It  I  - +  - 

-  14  - 


t 


(2.1a) 


and  =  c'<5£j»  where  d  is  the  layer  thickness,  and  c  (c.^,  ct)  and  p  (pd)  are 
the  acoustic (dilatational  d,  and  transverse  t)  wave  speeds  and  densities  in 
the  external  (or  layer)  media,  respectively.  The  quantities  C  and  C 

S  3 

depend  [6]  on  these  parameters  or  combinations  these  of,  such  as  on 


5  "  j  x  (ri\  -  sinzQ  ) 


(2.2b) 


where  x  «  27rf  d/c.  Viscosity  is  described  by  complex  wave  velocities  c^,  c  . 
For  the  case  of  a  fluid  layer  (  c  =  0),  one  has  C  =»  cot  5  ,  C  =>  tan  6  . 

,  t  S  cl 

A  taylor  expansion  of  the  denominators  of  Eqs.  (2.1)  leads  to  the  series 
of  resonance  terms  (of  Breit-Wigner  form) 


where 
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(2.4a) 


are  found  as  the  solutions  of  the 


The  (normalized)  resonance  frequencies  x 

s,a 

lossless  free-plate  characteristic  equations 


(2.5) 


and  the  resonance  widths  T  as  well  as  their  loss  broadening  AF  are 

ms,a  s,a 

given  [6,7]  as  analytic  expressions  in  terms  of  the  material  parameters  of 
the  plate.  As  mentioned,  this  leads  to  the  possibility  of  determining  these 
parameters  from  a  measurement  of  the  resonance  features. 

The  derivation  of  Eqs.  (2.3)  -  (2.5)  was  based  on  the  assumption  of  a 
large  impedance  contrast  between  fluid  and  plate,  i.e.,  T  «  1.  A  Mittag- 
Leffler  type  meromorphic  series  expansion  gives,  e.g.,  for  a  lossless  fluid 
layer  : 
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where  a  =  tanh  T, 
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Equation  (2.6)  is  an  exact  representation  for  T,  valid  for  all  T  <  1. 

A  corresponding  expression  can  be  found  for  T  >  !  [10]. 

3.  NUMERICAL  RESULTS. 

2 

Figure  3.1  shows  the  resonances  in  ] T J  vs.  f  at  8  =  19°,  for  a  one-inch 
thick  absorptive  plexiglas  plate  in  silicone  oil,  representing  plots  of  the 
individual  resonance  terms  in  the  present  formalism  [7].  The  classification 
in  terms  of  antisymmetric  (a)  and  symmetric  (s)  plate  waves,  as  well  as  their 
dilatational  (d)  and  transverse  (t)  content  follows  Brekhovskikh  [2]. 


Figure  3.2  then  shows  the  coherent  sum  of  the  terms  of  the  present  FMU  , 

resonance  theory,  compared  to  the  result  of  the  exact  expression,  Eq.  (2.1a). 

The  agreement  is  adequate  at  the  prevailing  impedance  contrast  (t  =  0.211  at  ' 
8  =  0°),  and  it  also  shows  that  the  representation  by  individual  resonances 
of  the  FMU  theory  provides  a  means  of  resolving  cases  of  interfering  resonances 
in  the  coherent  sum,  such  as  at  f  -  60,  105  and  155  kHz. 

4.  OCEAN  1R  RESONANCES.  - 

Cal:;. tated  reflection  coefficients  for  models  of  a  layered  ocean  floor 
often  contain  prominent  resonances  [15].  Figure  4.1  shows  a  two-dimensional 
plot  of  j  R |  for  two  (fluid)  sediment  layers  on  top  of  a  fluid  substratum, 
plotted  vs.  both  f  and  0  [14].  In  this  reference,  a  detailed  analysis  is 
given  how  to  determine  the  bottom  structure  from  measurements  on  the  resonances 
a  large  redundancy  of  measurements  was  found  to  apply.  In  addition,  the  ringing 
of  layer  resonances  as  manifested  by  a  tail  in  long,  specularly  reflected 
wavetrains  (Fig.  4.2)  has  been  suggested  [16]  as  a  means  of  measuring  the 
resonance  properties. 

5.  CONCLUSION. 

The  analysis  of  resonances  in  the  acoustic  reflection  from  layers  and 
plates  has  shown  that  reflection  and  transmission  coefficients  can  be 
represented  as  a  sum  of  resonance  terras,  of  Breit-Wigner  form  familiar  from 
nuclear  physics,  either  in  the  "resonance  approximation"  or  in  an  exact 
meromorphic  series.  The  resonance  positions  (in  frequency,  or  in  angle  of  . 

incidence),  heights,  and  widths  depend  on  the  material  parameters,  and  their 
measurement  may  thus  serve  to  determine  these  parameters  by  remote  acoustic 
sensing.  An  explicit  demonstration  of  such  a  procedure  was  provided  [14,16] 
for  the  case  of  sediment  layers  on  the  ocean  floor. 

Portions  of  the  results  reported  here  are  based  on  work  supported  by 
the  ONR  and  the  NSWC. 
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Fig.  3.1. 
Fig.  3.2. 
Fig .  4.1. 
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Individual  resonances  of  |t|  from  the  resonance  theory,  for  a 
1-inch  plexiglas  plate  in  silicone  oil  at  6  =  19°. 

2 

Coherent  sum  of  resonances  in  j  T |  for  a  I  -  inch  plexiglas  plate 
in  silicone  oil  at  0  =  19°,  compared  with  exact  result. 

2 

] R |  for  two  ocean  -  floor  sediment  layers  on  top  of  a  substratum, 
plotted  vs.  frequency  -  thickness  and  incident  angle. 


Fig.  4.2.  Reflected  signal  vs.  time,  for  incident  wave  train  with  carrier 

frequency  equalling  a  resonance  frequency  of  an  ocean  floor  sediment 
layer.  Final  transient  is  due  to  ringing  of  layer  resonance. 
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:  The  classical  normal-mode  series  of  acoustic  scattering  from  solid 

elastic  cylinders  and  spheres  is  reformulated  in  terms  of  the  S  -  function 
as  developed  in  nuclear  scattering  theory.  It  is  then  subjected  to  the 
Watson  transformation,  which  permits  an  evaluation  of  the  scattering 
amplitude  at  its  poles  ("Regge  poles")  and  saddle  points  in  the  complex 
mode-number  plane.  The  saddle  point  contributions  are  obtained  after 
expanding  the  amplitude  in  a  Debye  series,  and  correspond  to  a  reflected 
wave  and  to  transmitted  dilatational  and  shear  waves  that  undergo  internal 
reflections  and  mode  conversions.  The  theory  of  these  waves  was 
experimentally  verified  by  Quentin  et  al.  The  pole  residues  furnish 
circumferential  (surface,  creeping)  waves  which  are  of  both  Franz  type 
(propagating  externally),  and  of  elastic  type  (Rayleigh  and  Whispering 
Gallery  waves,  propagating  internally).  The  theory  of  these  waves  was 
experimentally  verified  by  Ripoche  et  al. 

:  On  reformule  la  serie  classique  des  modes  norraaux,  qui  decrit  la  diffusion 

acoustique  par  des  cylindres  et  des  spheres  solides  elastiaues  en  l'exprimant 
par  la  fonction  S, comme  elle  eat  developpee  dans  la  theorie  de  la  diffusion 
nucleaire.  Elle  est  alors  soumise  3  la  transformation  de  Watson,  ce  qui 
permet  une  evaluation  de  l'amplitude  de  diffusion  a  ses  poles  ("poles  de 
Regge")  et  3  ses  cols  dans  le  plan  complexe  du  numero  de  mode.  Les 
contributions  des  cols  sont  obtenues  apres  le  developpement  de  1' amplitude 
en  "serie  de  Debye",  et  elles  correspondent  3  une  onde  reflechie  et  3  des  ondes 
transmises  de  type  dilatation  et  cisaillement,  soumises  3  des  reflexions 
internes  et  des  conversions  de  mode.  La  theorie  de  ces  ondes  a  ete  verifiee 
experimentalement  par  Quentin  et  al.  Les  residus  des  poles  fournissent  des 


ondes  circonferentielles  (ondes  de  surface,  "Creeping  waves",  ou  du  type  de 
Franz,  se  propageant  sur  le  cote  externe  de  1 ' interface) ,  ou  du  type 
Slastique  (ondes  de  Rayleigh  et  de  galerie  A  echo, se  propageant  sur  le  cote 
interne  de  1 ’ interface)  .  La  theorie  de  ces  ondes  a  ete  verifiee 
experiment a lement  par  Ripoche  et  al. 

1.  INTRODUCTION. 

The  first  comprehensive  theory  of  sound  scattering  from  elastic 
cylinders  and  spheres  was  presented  by  Faran  [1]  in  1951  ;  further  extensive 
work  on  spheres  is  due  to  Hickling  [2].  In  their  analysis,  the  acoustic  field 
is  represented  as  a  Rayleigh  -  type  series  summed  over  normal  modes  labeled 
by  the  mode  number  n,  the  modes  having  been  obtained  by  the  separation  of 
variables  in  cylindrical  or  spherical  coordinates,  respectively.  Although 
this  solution  can  be  numerically  evaluated  to  obtain,  e.g.,  the  backscattering 
cross  section  as  a  function  of  frequency,  it  does  not  lend  itself  to  any 
physical  understanding  of  the  scattering  process. 

Such  an  understanding  was  provided,  however,  by  the  work  of  Franz  [3,4] 
on  the  scattering  of  electromagnetic  waves  from  perfect  conductors,  this 
being  the  analogue  of  acoustic  scattering  from  impenetrable  (soft  or  rigid) 
objects.  His  application  of  the  Watson  transformation  [5]  provided  a 
separation  of  the  scattered  field  into  two  parts  ,one  being  a  wave  which  in 
the  limit  of  large  values  of  ka  (with  k  the  acoustic  wave  number,  a  the 
radius  of  the  cylinder  or  sphere)  corresponds  to  geometrical  reflection 
from  the  apex  of  the  object,  and  the  other  representing  a  series  of  surface 
waves  (called  "creeping  waves"  by  Franz)  which  propagate  circumferentially 
around  the  scatterer,  thus  in  effect  making  up  the  diffraction  phenomenon. 

This  approach  was  carried  further,  and  applied  to  the  case  of  penetrable 
cylinders  [6-9]  and  spheres  [10,11].  Here,  interior  fields  are  present,  and 
the  following  separation  of  the  total  field  can  be  made  :  (a)  a  "geometric" 
part  which  in  the  large-ka  limit  corresponds  to  both  a  specular  reflection 
from  the  apex,  and  to  rays  refracted  into  the  interior  of  the  scatterer, 
which  then  re-emerge  into  the  exterior  fluid  either  immediately,  or  after 
a  series  of  multiple  internal  reflections.  For  an  elastic  object,  mode 
conversion  into  rays  of  shear  type  may  occur  during  these  refractions  or 
internal  reflections  ;  (b)  a  surface  wave  part  in  which  circumferential 
waves  propagate  around  the  scatterer  both  externally  (of  "Franz  type", 
similar  to  those  for  impenetrable  objects),  and  internally  (of  "Rayleigh" 
and  "Whispering  Gallery"  type,  the  former  corresponding  to  the  Rayleigh 
wave  in  the  limit  of  a  flat  surface). 
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Comprehensive  experimental  studies  have  been  performed  which 
demonstrated  qualitatively  and  quantitatively  the  correctness  of  this 
theory.  At  large  frequencies  (ka  „>  100)  where  the  geometrical  fields 
dominate,  studies  of  the  reflection  of  short  acoustic  pulses  from  elastic 
cylinders,  spheres  and  shells  were  carried  out  by  Quentin  et  al  [12  -  16]. 
They  verified  by  a  measurement  of  both  the  pulse  arrival  times  and  their 
amplitudes,  the  existence  and  the  predicted  properties  of  transmitted  waves, 
including  their  multiple  reflections  and  mode  conversions.  The  surface  waves 
on  elastic  cylinders  have  been  studied  extensively  by  Ripoche  et  al  [17  -  19] 
at  lower  values  of  ka  (<^50),  and  the  predicted  dispersion  and  absorption 
curves  of  both  Franz-type  and  elastic-type  (internal)  surface  waves  were 
confirmed  experimentally. 

2.  THEORY  OF  THE  SCATTERING  PROCESS. 

The  total  acoustic  field  in  the  presence  of  an  infinite  elastic 
cylinder  is  given  in  cylindrical  coordinates  (r,0,z)  by 

p=]T  (2-3^)^  ]  (kr)  \  cos-n.^  (2.1c) 

-n=o 


(corresponding  to  normal  incidence  of  sound  ;  for  oblique  incidence,  see  [20]) ; 
for  a  sphere,  it  is  in  spherical  coordinates  (r,9,0)  : 

oO 

P=  Y.  Urc+l)>t>v  {  ^(coiG).  (2. is) 

na=  O 

The  incident  pressure  field,  given  by  the  first  terms  in  Eqs.  (2.1),  is  here 
normalized  to  unity,  and  Tn  is  the  partial-wave  scattering  amplitude 
("T  -  function") in  the  normal-mode  or  Rayleigh  series  of  Eqs  (2.1).  It  is 
customary  in  nuclear  physics  to  rewrite  the  total  amplitude  in  the  form 

P  ”  jE(2-^u>)-i’'{H^,(U-r)+  <2-2c) 
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where 
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gives  the  relation  between  T  andthe"S  -  function"  S  . 
r  n  n 

Satisfying  the  boundary  conditions  at  r  =  a  leads  to 
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are  the  S  -  functions  for  a  soft  scatterer  (for  rigid  scatterers,  Hn  or 
hn^/  appear  ).  The  quantities  F  are  given  in  the  literature  11];  they  are 
proportional  to  p/PQ»  the  density  ratio  of  the  fluid  and  the  scatterer,  so 
that  F  «>  for  a  soft,  and  F  -*■  0  for  a  rigid  object. 

The  Watson  transformation  [3-5]  consists  in  rewriting  the  normal-mode 
sum  as  a  contour  integral  in  the  complex  n-plane  : 
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where  C  tightly  surrounds  the  positive  real  axis  in  the  complex  s>  or  X 
plane,  passing  through  v  ■  0  (P  »  principal  value)  but  to  the  right  of  X  *  0. 
The  Imai  transformation 
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where 


Q^(t>(c°s9)=-£  {P^ase)^  ~^Qr (2.9) 
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defines  Q  in  terms  of  the  Legendre  function  of  the  second  kind,  Q  , 

^  .  (2)  y 

splits  Eqs .  (2.7)  into  two  portions  (note  that  the  terms  with  H  or 

(2) 

1»X_  1/2  i-nte8rate  to  zero).  The  first  one,  containing  cos  v  ( 0  —  tt )  or 

wo  (cos  8),  can  be  evaluated  at  the  poles  of  the  S-function  S  in  the 

A-i/z  ( i )  v 

V  -  planeC'Regge  poles"),  given  by  the  zeros  of  -  z^  and  denoted  by 

v  =  (■£  =  1 ,2,3. .  .labeling  their  multiplicity)  : 
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where  S  is  the  residue  of  S  at  V  =  The  exponential  form  of 
cos  V  -  (0  -  tt)  ,  and  the  asymptotic  form 


U.llj 


demonstrate  that  Eqs.  (2-10)  represent  circumferential  waves  pcw  with 
propagation  constant  for  the  cylinder,  +  1/2  for  the  sphere.  From 
this,  one  finds  the  phase  velocities 
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and  the  attenuation  angles  (amplitude  exp 
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of  the  surface  waves,  both  the  external  and  the  internal  ones. 
The  second  part  of  Eqs.  (2.7)  is 
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These  integrals  can  be  evaluated  at  the  saddle  points  of  the  integrand 

[3, 4, 8, 9]  and  are  then  seen  to  represent  geometrical  waves  p  ,  either 

gw 

reflected  from  the  apex  of  the  scatterer  or  undergoing  internal  transmissions 
without  and  with  multiple  internal  reflections,  including  mode  conversions 
for  an  elastic  scatterer.  This  is  obtained  from  an  expansion  of  in 
Eqs.  (2.13)  into  a  Debye  series,  as  was  done  for  an  elastic  cylinder  by 

ii 

Brill  and  Uberall  [9],  but  which  for  the  simpler  case  of  fluid  scatterers 
[8,10]  becomes,  e.g.  for  cylinders  : 


^  fS)(V  _  ^ 

^  l^±T. 
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where  6  =  c/cq,  the  sound  velocity  ratio  of  fluid  and  scatterer.  This 
contains  an  external  reflection  coefficient 
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and  an  internal  one. 


(2.14b) 


as  well  as  the  transmission  coefficients 
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(2.14c) 
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into  and  out  of  the  object,  respectively  ;  also,  N  =  Sp/p  •  For  the  sphere, 
one  replaces  1  -*•  h^  1  in  Eqs.  (2.14).  The  individual  terms  in  Eqs.  (2.14), 
evaluated  at  their  saddle  points  in  Eqs.  (2.13),  furnish  the  geometrical  rays 
corresponding  to  no  (R^)  or  ^  internal  traversals. 

3  EXPERIMENTAL  RESULTS. 

The  surface  waves  form  resonating  standing  waves  around  the  scatterer  for 
Re  V-  -  n,  where  n  (n+1/2)  of  their  wavelengths  span  the  cylinder  (sphere  ; 
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here,  a  X/4  phase  jump  cakes  place  at  each  of  their  two  convereence  points  [21]). 
By  observing  these  resonances  on  cylinders,  Ripoche  et  al  [17-19]  have  verified 
the  surface  waves,  and  determined  their  phase  velocities  c£  and  attenuations 
Eqs.  (2.12).  Figure  3.1  shows  their  results  (crosses)  for  the  attenuation 
of  Whispering  Gallery  waves  (internal,  ,£  =  2,3)  as  compared  to  theory  [7] 

(curves) . 

For  the  transmitted  waves  in  cylinders  and  spheres,  experiments  [12,13] 
give  excellent  agreement  with  the  theoretical  arrival  times  and  amplitudes  [9]. 
Figure  3.2  shows  measured  and  calculated  arrival  times  of  short  acoustic  pulses 
backscattered  from  a  lucite  cylinder  (the  inserts  show  the  type  of  traversals). 

4.  CONCLUSION 

The  theory  of  sound  scattering  from  elastic  cylinders  and  spheres,  as 
based  on  the  Watson  transformation,  gives  an  accurate  account  of  the  physics  of  the 
diffraction  phenomenon.  Its  various  aspects,  i.e.  the  existence  and  properties 
of  (resonating)  external  and  internal  circumferential  waves,  and  of 
geometrically  reflected  and  internally  transmitted  rays,  have  all  been  verified 
experimentally.  Analogous  studies  with  scatterers  of  more  general  shapes  are 
now  called  for. 

Portions  of  the  results  reported  here  are  based  on  work  supported  by  the 
Office  of  Naval  Research. 
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FIGURE  CAPTIONS. 

Fig.  3.1.  Attenuation  of  the  whispering  Gallery  waves  £  =  2,  3. 

Solid  curves  :  theory  [7]  ;  crosses  :  experiment  [19]. 

Fig.  3.2.  Experimental  [12]  and  theoretical  [9]  arrival  times  of  short 
sound  pulses  traversing  a  lucite  cylinder  (Solid  rays  : 
dilatational  broken  rays  :  shear). 
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The  eigenf requencies  of  a  cylindrical  obstacle  (of  finite  or  infinite 
length)  can  be  interpreted  as  the  resonances  due  to  phase  matching  of 
circumnavigating  helical  surface  waves  .  For  the  case  of  a  cylinder  of  finite 
length,  the  pitch  angle  of  the  helix  can  assume  a  discrete  set  of  values 
only.  Resonant  eigenvibrations  can  be  excited  by  waves  incident  in  an  oblique 
fashion,  which  generates  the  helical  waves.  A  refraction  effect  is  found  to 
take  place  between  the  incident  and  the  helical-wave  directions.  We  obtain 
pole  diagrams  of  the  scattering  amplitude  in  the  complex-frequency  plane, 
by  using  the  T-matrix  approximation  for  finite  cylinders.  In  addition,  pole 
diagrams  for  spheroidal  scatterers  are  obtained  by  the  use  of  the  T-matrix 
and  of  spheroidal  wave  functions.  While  the  poles  of  symmetric  scatterers 
(spheres  or  infinite  cylinders)  are  degenerate  in  the  azimuthal  quantum 
number  m,  the  degeneracy  for  the  poles  of  finite  cylinders  and  of  spheroids 
is  lifted.  This  m-splitting  is  explained  by  the  phase  matching  of  helical 
waves  with  various  allowed  pitch  angles.  Dispersion  curves  for  the  phase  and 
group  velocities  and  attenuations  of  the  helical  waves  are  obtained. 

Les  frequences  propres  d'un  obstacle  cylindrique  (de  longueur  finie  ou 
infinie)  peuvent  etre  interpretdes  cornme  dues  a  1 1  accord  entre  les  phases 
d'ondes  se  propageant  sur  la  surface  d'une  fa<;on  helicoidale.  Dans  le  cas 
d'un  cylindre  de  longueur  finie,  l' angle  de  pas  de  1 'helice  ne  peut  prendre 
qu'une  serie  de  valeurs  discretes.  Des  vibrations  propres  resonnantes  peuvent 
Itre  excitees  par  des  ondes  incidentes  de  direction  oblique,  ce  qui  produit 
les  ondes  helicoidales .  Un  effet  de  refraction  est  trouve  entre  les  directions 
de  l'onde  incidente  et  de  l'onde  helicoi’dale.  On  obtient  des  diagrammes  de 
p81es  dans  l'araplitude  de  diffusion  dans  le  plan  complexe  de  la  frequence,  par 
un  calcul  ut i  1  isant  1.' aonroximation  de  la  matrice  T  pour  des  cylindres  <:inis. 
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En  plus,  on  obtient  des  diagramraes  de  poles  pour  des  obstacles  spheroi'daux 
en  utilisant  la  matri.ce  T,  ou  des  fonctions  d'ondes  spheroi'dales  .  Tandis  que 
les  p51es  d'obstacles  syraStriques  (sph&res,  ou  cylindres  inf inis )  degenerent 
vis-a-vis  du  nombre  quantique  azimuthal  m,  cela  n'est  plus  le  cas  pour 

les  poles  de  cylindres  finis  et  de  spheroides .  La  separation  rdsultante  entre 
les  valeurs  de  m  s'explique  alors  par  l'accord  de  phases  des  ondes 
heli  coi'dales  possedant  diffdrents  angles  d' inclinaison  permis.  On  obtient 
des  courbes  de  dispersion  pour  les  vitesses  de  phase  et  de  groupe  des 
ondes  hdlicoi'dales  . 

1.  INTRODUCTION. 

The  most  general  case  of  the  scattering  of  waves  obliquely  incident  on 
an  infinite  cylinder  has  first  been  considered  by  White  [1]  ;  he  deals  with 
a  plane  wave  in  an  elastic  medium  incident  on  an  elastic  cylindrical 
inclusion.  This  theory  has  recently  been  completed  by  Delsanto  et  al  [2,3]  ; 
it  is  based  on  the  classical  normal-mode  (or  partial-wave)  Rayleigh  series 
expansion.  Special  cases  have  been  investigated  in  greater  detail  :  acoustic 
waves  incident  on  an  elastic  cylinder  [4]  or  shell  [5],  or  on  a  fluid  cylinder 
[6,7],  and  elastic  waves  incident  on  a  fluid-filled  cylindrical  cavity  [8] 

(in  a  similar  fashion,  elastic  waves  on  a  spherical  cavity  have  also  been 
studied  [9]).  In  these  latter  investigations,  the  complex  eigenfrequencies 
of  the  mentioned  obstacles  have  been  explicitly  calculated,  as  well  as  the 
excitation  of  the  corresponding  eigenvibrations .  While  analogous  calculations 
cannot  be  carried  out  for  the  case  of  cylinders  of  finite  length,  due  to  the 
non-separability  of  this  problem,  results  for  the  eigenfrequencies  can 
nevertheless  be  obtained  here  by  the  use  of  special  methods .  If  only  the 
interior  vibrations  of  a  fluid  cylinder  in  vacuo  are  considered,  the  problem 
is  still  exactly  soluble,  and  the  corresponding  eigenvibrations  have,  for  the 
first  time,  been  interpreted  as  being  due  to  the  phase  matching  of  helical 
surface  waves  [10]  (of  internal  type,  in  this  case).  For  the  non-separable 
exterior  problem  of  finite-length  impenetrable  cylinders,  the  T-raatrix 
method  of  Waterman  [11]  has  been  employed,  modified  so  as  to  furnish  complex 
eigenfrequencies  [12].  A  very  comprehensive  study  of  this  problem  has  recently 
been  carried  out,  on  the  basis  of  this  method,  which  is  contained  in  a  paper 
[13]  that  also  di  scusses  the  acoustic  eigenfrequencies  of  impenetrable 
spheroids  obtained  by  the  use  of  spheroidal  functions.  In  this  study,  the 
phase  matching  of  external  helical  surface  waves  has  been  invoked  as  an 
explanation  for  the  finite-cylinder  eigenfrequencies,  and  for  their  splitting 


into  components  corresponding  to  different  values  of  the  azimuthal  quantum 
number  m  (the  latter  being  a  measure  for  the  pitch  angles  of  the  helical  waves 
which  due  to  the  finite  cylinder  length  form  a  discrete  set).  As  to  the 
excitation  of  helical  surface  waves  by  incident  sound,  if  was  found  [6,13] 
that  this  takes  place  in  a  refractive  way,  the  helical  pitch  angle  being 
different  from  the  incident  angle.  Experiments  on  the  scattering  of  obliquely 
incident  sound  by  elastic  cylinders  have  recently  been  initiated  [14]. 


INTERNAL  HELICAL  WAVES. 


A  fluid  cylinder  in  vacuo,  of  radius  a  and  length  Ladmitsan  internal 
acoustic  field  (wave  vector  k,  with  k  =40/ c)  : 


p  (^)  =  ( jcr)  e 


(2.1) 


where 


~K‘ 


-  kz  -  u:  ; 


(2.2a) 


the  single  -  valuedness  condition  under  (J>  -►  <p  +  2 ir  gives 

kf=  ,  (2.2b) 

and  the  boundary  condition  at  the  end  faces  leads  to 

!  I—  ^  *  4- >  ^  >  •  ••  •  (2.2c) 

The  Dirichlet  boundary  condition  at  r  =  a  gives 

n=i,2}...  (2. 2d) 

where  x.  is  the  nth  zero  of  J  (x) .  [The  Neumann  condition  for  a  fluid 
mn  m 

cylinder  in  a  rigid  enclosure  would  lead  to  the  zero  of  J'm(x)].  Inserted  in 

Eq.  (2.2a),  this  gives  the  eigenvalues  k  .  of  k. 

.  .  nmj 

We  now  introduce  a  tangential  wave  vector  k  of  the  surface  field,  with 


which  describes  the  propagation  of  helical  surface  waves.  The  conditions 

(2.2b)  and  (2.2c)  then  represent  the  phase  matching  of  such  waves  after 

circumnavigating  the  cylinder,  and/or  getting  reflected  from  the  end  faces. 

Equations  (2.2b),  (2.2c)  inserted  in  Eq.  (2.3)  furnish  eigenvalues  (k  )  . 

s  nij  • 

The  helical-wave  phase  velocities  c  =  oj/k  can  be  obtained  from 

s  s 

C^/c  m  k / kg  at  the  discrete  points  (resonance  frequencies)  where  phase 
matching  is  satisfied  :  ' 

</r)l(.a/l-)2]/(V+  tfr)x(xi/ L)1]}  (2.4) 

Note  that  these  correspond  to  helical  waves  of  pitch  angle  a  with  the 
z  -  axis,  tan  a  =  k^/kz,  which  due  to  the  finite  length  of  the  cylinder  assumes 
the  discrete  values 


tan  <K«i  -  (W//-T)  (L/-a). 


A  given  helical  wave  thus  corresponds  to  a  fixed  ratio  m/j.  For  an  infinite 
cylinder  (L  -*■«>,  j  -*•«)  a  is  continuous. 

Figure  2.1  shows  how  the  discrete  points  of  Eq .  (2.4),  corresponding  to 
the  eigenfrequencies  of  the  cylinder,  when  connected  according  to  Eq.  (2.5) 
furnish  the  dispersion  curves  of  the  helical  surface  waves. 

3.  REFRACTION  EFFECT. 

For  a  cylinder  in  a  medium,  helical  waves  can  be  generated  by  an 
incident  plane  acoustic  wave.  If  the  latter  arrives  at  an  angle  y  with  the  z- 
axis  so  that  k  =  k  siny,  k  =  k  cosy,  the  total  field  for  an  infinite 

y  z 

cylinder,  given  by  Eq.  (2.2c)  of  Ref.  [15],  gets  modified  to 


leading  after  application  of  the  Watson  transformation  to  the  creeping-wave 
sum  analogous  to  Eq.  (2.10c)  of  [15]  : 


35  - 


pew 


=  7re*  - 7T  ^  (kyr)  cosy  (#-*).  (3.2) 


‘ —  si-nTrV. 

J=1  1 


The  phase  factor  =»  z  +  <J>  Re  shows  that  these  surface  waves  are 

helical,  with  wave  fronts  a<J>  =  -(ak  /Re  V*  )  z  whose  normals  make  an  angle 
-1  z  X> 

$n  ~  tan  (Re  V^/ak^)  with  the  z-  axis.  This  defines  the  law  of  refraction 


ian  fa  =  la ny9  ^  -  Re 

between  incident  direction  y  and  helical-wave  direction  §p.  The  phase 
velocity  v^  =  ck/fcg  of  the  helical  waves  is 

U  pk  =  c/{  (He  tf/kxa-)1  f  cos2-^  }  J 


for  the  case  of  external  waves  on  rigid  or  soft  cylinders  where  the 
asymptotic  expansion  of  Franz  [16]  for  v^(k^a)  can  be  used,  one  has 


(3.3) 


(3.4) 


vk  *  - 


Sin^y  \2A 
ltd  ) 


(3.5) 


IT  S 

where  =  1.469354,  q^  =  3,372134.  Including  higher  terms,  Fig.  3.1  shows 
vjp  of  helical  waves  aty=0®  and  45°,  and  the  refraction  angle  <p^  for  y  =  45® 
for  a  soft  cylinder. 

4.  COMPLEX  EIGENFREQUENCIES . 

The  eigenfrequencies  for  a  cylinder  in  vacuo,  Section  3,  are  real 

since  no  radiation  loss  can  occur.  We  present  figures  showing  examples  of 
complex  eigenfrequencies  that  correspond  to  the  resonances,  due  to  phase 
matching,  of  external  circumferential  waves  on  elongated  objects. 

Figure  4.1  shows  the  eigenfrequencies  in  the  complex  kb  plane  of  an 
acoustically  rigid  prolate  spheroid,  for  various  axis  ratios  b/a  as  indicated 
(a  being  the  serairainor  axis).  These  were  obtained  [13]  by  satisfying  the 


boundary  condition  with  spheroidal  wave  functions.  One  notices  that  the 
sphere  values  (crossed  squares),  degenerate  in  m,  split  into  branches  non¬ 
degenerate  in  m  for  spheroids.  Figure  4.2  shows  the  same  for  the  magnetic 
eigenvibrations  of  a  conducting  cylinder  of  length  L  (the  electromagnetic 
analogue  of  an  acoustically  rigid  cylinder)obtained  with  the  T-matrix 
method  [11-13],  A  phase-matching  model  for  surface  waves  has  been  developed 
here  [13]  in  order  to  show  that  the  m-splitting  of  the  eigenfrequencies 
corresponds  to  helices  of  different  pitch  angles. 

5.  CAVITIES. 

Analogous  results  were  obtained  for  infinite  cylindrical  cavities,  using 
our  general  theory  [2,3].  Figure  5.1  shows  the  eigenfrequencies  in  the  complex 
k^a  plane  for  an  empty  cavity  in  aluminum,  for  both  compressional  (p)  and 
shear  (s)  type  surface  waves  (here,  k  is  the  dilatational  propagation  constant). 

6.  CONCLUSION. 

The  complex  eigenfrequencies  of  f ini te- length  cylinders  show  a  splitting 
according  to  the  azimuthal  quantum  number  m.  They  can  be  interpreted  as  the 
resonances,  due  to  phase  matching,  of  helical  surface  waves  of  different  pitch 
angles.  For  infinite  cylinders,  a  continuum  of  pitch  angles  occurs.  The  helical 
surface  waves  can  be  excited  by  incident  acoustic  waves,  and  refraction  takes 
place  between  the  incident  and  the  surface  wave  directions.  Recent  experiments 
[14]  are  now  investigating  these  problems.  Theoretically,  there  have  also  been 
geometrical  investigations  of  helical  waves  on  cylinders  [17], 
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FIGURE  CAPTIONS. 


Fig.  2. 
Fig.  3. 

Fig.  4. 
Fig.  4. 
Fig.  5. 


1.  Dispersion  curves  of  internal  helical  surface  waves  on  a  finite 
fluid  cylinder  in  vacuo,  of  dimension  1:1  (L  =  2a). 

1.  Dispersion  curves  of  external  helical  surface  waves  on  an  infinite 
soft  cylinder  for  y  =  90°  (dashed)  and  45°  (solid),  and  refraction 

angle  at  y  =  45°  (right  scale),  plotted  vs.  k  a. 

y 

1.  Complex  eigenfrequencies  for  a  rigid  spheroid  with  semi-major  axis  b. 

2.  Complex  magnetic  eigenfrequencies  of  a  conducting  cylinder  of  length  L. 

1.  Eigenfrequencies  in  the  complex  k  a  plane  for  an  empty  cylindrical 
cavity  in  aluminum.  ^ 
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SUMMARY  :  We  consider  the  scattering  of  elastic  dilatational  and  shear  waves  from 

cylindrical  and  spherical  cavities  and  inclusions  in  an  elastic  medium.  The 
normal  mode  series  of  the  scattering  amplitude  is  reformulated  in  terms  of 
the  S-  function,  and  the  poles  of  the  S-  function  in  the  complex  frequency 
plane  are  identified.  The  amplitude  is  rewritten  as  a  "background  term" 
including  specular  reflections  and  external  surface  waves,  plus  a  series  of 
(internal)  resonance  terms.  This  formulation  is  termed  the  "Resonance 
Scattering  Theory”  (RST) .  The  connection  between  the  resonances  and  the 
surface  waves  is  established  via  expressing  the  complex-frequency  poles  of 
the  scattering  amplitude  by  the  Regge  poles  in  the  complex-mode  number  plane, 
and  the  frequency  resonances  in  successive  modes  are  recognized  as  the  Regge 
recurrences  of  surface-wave  resonances.  This  permits  us  to  obtain  the 
dispersion  curves  of  phase  and  group  velocities  of  the  (internal)  surface 
waves  from  the  eigenf requencies  of  the  cavity.  We  also  mention  experiments 
on  ultrasonic  scattering  from  cavities  and  inclusions. 

RESUME  :  Nous  etudions  la  diffusion  d'ondes  elastiques  de  type  dilatation 

ou  cisaillement  par  des  cavites  et  inclusions  cylindriques  ou  spheriques. 

La  serie  des  modes  norrnau:;  donnant  l'amplitude  de  diffusion  est  recoraposee 
en  utilisant  la  fonction  S,  et  les  poles  de  la  fonction  S  dans  le  plan 
complexe  des  frequences  sont  identifies.  L '  ampl  i  tude  es  t  reecr  i  te  sous  la  forme  d  'un 
"terme  de  fond",  comprenant  reflexions  speculaires  et  ondes  de  surface 
externes,  plus  une  serie  de  termes  de  resonances  (internes).  Ce  formalisme 
est  appele  "  Theorie  de  diffusion  resonnante" ["Resonance  Scattering  Theory" 
(RST)].  Le  rapport  entre  les  resonances  et  les  ondes  de  surface  est  etabli 
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en  exprimant  les  poles  en  frequence  complexe  de  1 'amplitude  de  diffusion 
par  les  poles  de  Regge  dans  le  plan  complexe  du  numero  de  mode,  et  les 
resonances  en  frequence  dans  les  modes  successifs  sont  reconnues  come  les 
recurrences  des  resonances  d'ondes  de  surface.  Cela  nous  permet  d'obtenir 
les  courbes  de  dispersion  pour  la  vitesse  de  phase  et  de  groupe  des  ondes 
(internes)  de  surface,  en  utilisant  les  frequences  propres  de  la  cavite. 

Nous  mentionnons  egalement  des  experiences  sur  la  diffusion  ultrasor.ore 
par  des  cavites  et  inclusions. 

1.  INTRODUCTION 

The  Resonance  Scattering  Theory  (RST)  has  first  been  formulated  ([1]  ; 
see  also  [2,3])  for  the  case  of  sound  scattering  from  a  solid  cylinder 
or  sphere  in  a  fluid.  The  resonance  frequencies,  which  may  e.g.  be  represented 
by  a  "pole  diagram"  in  the  complex  frequency  plane  [4,5],  were  shown  to 
originate  from  the  phase  matching  of  surface  waves  [2,3,6].  When  the 
frequency  poles  are  transferred  to  the  complex  mode  number  plane  [6],  where 
they  constitute  the  "Regge  poles",  they  determine  by  their  real  and  imaginary 
parts  the  propagation  constant  and  attenuation  of  the  surface  waves, 
respectively  ;  hence,  the  dispersion  and  attenuation  curves  of  the  surface 
waves  can  be  obtained  from  the  resonance  frequencies  without  having  to  go 
through  the  (complex)  Watson  transformation. 

A  similar  program  may  be  carried  out  for  ultrasonic  wave  scattering  from 
cavities  of  cylindrical  [7-9]  or  spherical  [10-13]  shape,  or  from  solid 
inclusions  [14,15].  This  topic  is  chosen  here  to  present  the  resonance  theory, 
since  surface  waves  and  resonances  in  acoustic  scattering  from  elastic  objects 
are  treated  elsewhere  in  this  issue  [16,17].  While  experiments  on  resonant 
acoustic  scattering  from  solids  are  quite  numerous  by  now  (see,  e.g.,  [18,19] 
and  references  quoted  therein),  ultrasonic  experiments  on  scattering  from 
cavities  or  inclusions  as  done  so  far  [20-23]  have  not  been  specifically 
designed  to  detect  resonances,  so  that  resonance  effects  appear  in  their 
results  only  indirectly. 

2.  RESONANCE  THEORY. 

We  consider  an  infinite-cylindrical  cavity  in  a  solid  (density  p,  wave 
speeds  c  ,  c  )  filled  with  a  fluid  (density  p  .sound  speed  Cf)-  The 
displacement  field  u  in  the  solid  is  represented  by  a  scalar  (<J)  and  a 
vector  (4*)  potential  : 

AX.  =  V$  +  Vx  Y  (2.i) 
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(with  V.  ¥  =  0)  ,  and  a  normally  incident  p-wave  is  given  by 


=£  (2-S  U-Jn(J^)c«**, 


(2.2) 


I  vv  c- 
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where  k ^  =  tu/c^  (suppressing  exp-itot) .  Scattered  elastic  waves  are  of  both 


p-type. 
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or  s-  type  ("mode  conversion"). 
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where  k  =  co/ c  .  A  standing  wave  exists  in  the  fluid  filler  : 
s  s 
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where  =  o/c^.  Similar  equations  hold  [8]  for  s-wave  incidence  (assumed  SV, 

so  that  only  f  ;  f  ^  0) ,  with  coefficients  S  Sp,  S  SS  and  C  S.  The  sets  of 
J  z  '  n  n  n 

S  form  the  2x2  S-matrix. 
n 


Satisfying  the  boundary  conditions  (continuity  of  normal  stresses  and 
displacements,  zero  tangential  stress)  at  the  cylinder  radius  r  =  a  leads 
to  3x3  linear  systems 
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where  D  ^’S  are  3x3  matrices  and  E  ^,S  vectors,  all  with  known  elements  [8 
n 

which  consist  of  cylinder 
B  =  k  a  and  y  =  k  a,  and 


which  consist  of  cylinder  functions  of  order  n  and  arguments  a  =  k^a. 


Using  Cramer’s  rule*  the  elements  S  are  thus  fractions  with  D^Cto)  -  det(Dn  ) 
s 

=  -  det  (D  )  in  the  denominator.  Resonances  occur  when  D  (to)  =  0  ;  the 
n  '  n 

complex  solutions  to  .  (j  =  1,2,...)  represent  resonance  frequencies  common  to 

^  ilc  .  • 

all  four  amplitudes  .  Alternately,  at  a  given  to  (real),  (to)  =  0  is  an 

equation  for  n  with  complex  solutions  n  =  V  ^  (to) .  While  to^  are  complex- 

frequency  poles  ("SEM  poles",  see  [17])  of  the  scattering  amplitude,  v.  (to)  are 

its  complex  mode  number  ("Regge")  poles  [16]  ;  both  are  related  through 

D  (to)  =0. 

•  •  •  ik 

More  explicitly,  one  may  write  for  ,  e.g.. 
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with  Z  ^  expressed  by  the  elements  of  D  **,  and  S  -  eyp  2 i 6  the  value 

n  n  n  ^  n 

of  S/P  for  an  empty  cavity  (p^  =  0).  Expanding  (a)  about  the  resonance 

value  a  .  =  to  .a/c  leads  to 
nj  nj  p 
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This  is  the  Breit-Wigner  form  of  the  scattering  amplitude,  as  developed  in 
nuclear  scattering  [24].  It  consists  of  a  sum  of  resonant  terms  (with  ^nj* 
internal  resonances) ,  and  a  smooth  background  term  corresponding  to  the 
scattering  from  an  empty  cavity,  second  term  in  Eq.  (2.8).  However,  this 
background  still  contains  resonances  corresponding  to  external  (Franz)  waves, 
as  well  as  geometrical  (specurlarly  reflected,  and  transmitted)  waves  [16]. 

The  quantities  v^.  (w)  are  the  propagation  constants  of  surface  waves, 
as  in  the  Watson-transform  method  [16].  This  can  be  seen,  e.g.,  by  expressing 
the  resonance  denominator  in  CX  in  Eq.  (2.8)  by  one  linear  in  n  (assumed  a 
continuous  variable)  [6,  25]  ;  this  furnishes  the  Watson  results  without 
having  to  go  through  the  complex  arithmetic.  The  surface  wave  parameters  are 
thus 


sCj  ^  =  xi to/ Re.  Vj 

(phase  velocity) , 

sQ.  (oj)  =  -a/ Re  (dvj/dui) 

(group  velocity),  and  (amplitude  exp  -  : 

ft  =  i/lm ^ 

(attenuation).  Their  wavelength  (in  the  fluid)  is 


(2.10a) 
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showing  that  a  resonance  (Revj  ■+  n)  occurs  when  the  circumference  of  the 
cavity  equals  an  integral  number  of  wavelengths  of  the  surface  wave  ("phase 
matching") . 
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For  the  internal  surface  waves  which  have  small  imaginary  parts,  one 
has  at  resonance  (to  =  Reojnj)  approximately  : 


/CjP  (  Ra  ^ j  )  = 

/CjW  (Re  =  -cl  dPe  cu^j  /dn 


(2.11a) 


(2. 1  lb) 


(to  be  understood  as  a  difference  quotient),  so  that  the  dispersion  curves 
of  surface  waves  at  the  resonance  points  can  be  directly  read  off  trom  the 
latter . 

A  similar  theory  was  developed  for  the  case  of  scattering  from  spherical 
[14]  and  cylindrical  [15]  inclusions. 

3.  NUMERICAL  RESULTS. 

The  (internal)  frequency  resonances  appear  in  a  plot,  e.g.  of  |  SnPP  -  l|, 
vs.  k^a,  as  shown  in  Fig.  3.1(a)  for  a  water-filled  cylindrical  cavity  in 
aluminum.  Here,  the  resonance  and  background  terms  of  Eq.  (2.8)  interfere  ; 
but  if  the  soft-cylinder  background  is  coherently  subtracted  from  S^PP,  one 
obtains  the  pure  resonances  of  Fig.  3.1  (b) .  One  notices  that  peaks  of  the 
same  order  j  recur  in  successive  modal  amplitudes  at  successively  higher 
frequencies  ;  they  represent  "Regge  recurrences"  of  the  jth  surface  wave 
spanning  the  cavity's  circumference  successively  with  n  =  1 ,2 , 3 .. .wavelengths . 
The  corresponding  phase  velocity  dispersion  curves,  Eq.  (2.10a),  are  shown  in 
Fig.  3.2. 

The  background  terms  S  ^  ^  for  an  empty  cavity  themselves  have 

/  \ 

resonances,  but  with  large  imaginary  parts  so  that  e.g.  curves  I °  ^  ~  H 
appear  smooth,  as  the  broad  lobes  underlying  the  curves  of  Fig.  3.1(a).  One 
may  find  their  frequency  resonances,  however,  by  a  complex  pole  search,  as 
shown  in  Fig.  3.3  for  a  spherical  void  in  steel  [13].  The  poles  in  the 
a  =  k^a  plane  clearly  fall  into  layers  due  to  Rayleigh  (R),  compressionnal  (P) , 
and  shear  (S)  type  (external)  surface  waves  ;  their  dispersion  curves  are 
shown  in  Fig.  3.4.  Comparing  this  with  Fig.  3.2,  one  sees  that  dispersion 
curves  for  external  surface  waves  start  from  the  origin  and  reach  a  plateau 
for  high  frequencies,  while  for  internal  surface  waves,  they  start  at  infinity 
from  a  low-frequency  cutoff,  whence  they  descend  to  a  plateau.  This  can  be 
understood  geometrically. 

For  a  solid  inclusion  in  a  solid,  results  are  shown  in  Figs.  3.5  and  3.6. 
If  the  total  scattering  amplitude,  e.g.  the  modulus  of  Eq.  (2.3),  is  plotted 
vs.  frequency,  a  rather  smooth  curve  results  (solid  curve  in  Fig.  3.5  [22]). 


If,  however,  one  plots  individual  modes  n  in  the  fashion  of  Fig.  3J  (b) , 
peaks  appear  (depending,  of  course,  on  the  materials),  which  may  be  used  to 
obtain  dispersion  curves,  Fig.  3.6.  These  indicate  the  existence  of  two  types 
of  surface  waves  (solid  and  dashed) ,  probably  also  a  Rayleigh  wave  (lowest 
dashed  curve) . 

4.  EXPERIMENTS. 

Experimental  results  on  ultrasonic  wave  scattering  from  cavities  and 
inclusions  were  obtained  both  in  the  frequency  domain  and  in  the  time  domain 
(short  pulses)  [20  -  23].  Pulse  experiments  obtained  geometrical  and  surface 
wave  paths  from  arrival  times.  Besides  spherical  and  cylindrical  shapes, 
spheroids  have  also  been  investigated.  The  broken  curve  in  Fig.  3.5  shows 
backscattering  results  for  a  spherical  tungsten  carbide  inclusion  in 
titanium  [22],  as  compared  to  theory. 

5.  DISCUSSION. 

The  resonance  scattering  theory  (RST)  as  developed  for  acoustic  scattering 
from  solids,  has  here  been  applied  to  ultrasonic  scattering  from  cavities  and 
inclusions.  Internal  resonances  appear  prominently  if  the  properties  (e.g. 
density)  of  cavity  material  and  its  filler  differ  substantially  ;  if  this  is 
not  the  case,  they  are  less  evident.  External  resonances  on  cavities  are 
always  hard  to  see  ;  but  in  all  cases  (internal  and  external),  resonance 
frequencies  in  the  complex  plane  can  be  calculated  no  matter  how  large  their 
imaginary  part,  and  their  layers  define  clearly  the  various  types  of  surface 
waves.  Experiments  have  not  yet  been  designed  specifically  to  detect  resonances; 
in  a  sense,  the  pulse  arrivals  indicate  these  indirectly,  being  a  coherent 
sum  of  a  number  of  resonances  [17]. 

Portions  of  this  work  were  supported  by  the  Office  of  Naval  Research,  and 
by  the  Army  Research  Office. 
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FIGURE  CAPTIONS 

Fig.  3.1.  Modal  resonances  for  a  water-filled  cylindrical  cavity  in 
aluminum  :  (a)  jSnPP-l|  plotted  vs.  k^a, 

(b)  | SnPP  -  Sn(°)pp-lj,  background  subtracted. 

Fig.  3.2.  Dispersion  curves  of  internal  surface  waves  in  a  water-filled 
cylindrical  cavity  in  aluminum. 

Fig.  3.3.  Complex-frequency  poles  of  scattering  amplitudes  from  a  spherical 
void  in  steel,  showing  resonances  due  to  Rayleigh  (R), 
compressional  (P)  and  shear  (S)  type  surface  waves. 

Fig.  3.4.  Dispersion  curves  of  surface  waves  whose  resonances  furnish 
the  poles  in  Fig.  3.3. 

Fig.  3.5.  Theoretical  (solid)  and  experimental  (broken  curve)backscattering 
from  a  spherical  tungsten  carbide  inclusion  in  titanium. 

Fig.  3.6.  Dispersion  curves  for  two  families  of  surface  waves  on  an  iron 
sphere  in  aluminum. 
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SUMMARY  :  A  microinhomogeneous  medium,  consisting  of  randomly  distributed  cavities 
or  inclusions  in  a  homogeneous  elastic  matrix,  can  be  represented  as  a 
dispersive  homogeneous  medium  with  effective  material  constants  (moduli, 
bulk  wave  speeds,  and  absorptions).  For  wavelengths  long  compared  to  the 
size  of  the  scatterers,  Kuster  and  Toksoz  have  developed  a  method  (not 
including  rescattering)  which  obtains  these  effective  material  properties 
by  comparing  exact  and  effective  monopole,  dipole  and  quadrupole  amplitudes. 
We  extend  this  approach  to  the  case  where  the  wavelength  is  comparable  to 
the  size  of  the  scatterers  (assumed  spherical)  ;  in  this  case,  particle 
resonances  are  taken  into  account  and  lead  to  widened  resonances  in  the 
effective  material  parameters.  The  cases  of  bubbly  liquids,  of  perforated 
solids,  and  of  solids  with  solid  inclusions  (particulate  composites)  are 
treated  in  this  fashion.  Measurements  by  Kinra  and  Anand  verify  our  results. 
In  addition,  many  previous  results  for  the  effective  moduli  of  composites, 
obtained  in  the  static  (i.e.,  low-frequency)  limit,  are  recovered  as 
particular  cases  of  our  approach. 

RESUME  :  Un  milieu  micro-inhomogene,  qui  consiste  en  cavites  ou  inclusions 

distributes  dans  une  matrice  elastique  homogene  de  fa<;on  aleatoire, 
peut  etre  considere  comme  un  milieu  homogene  et  dispersif  possedant  des 
constantes  de  materiaux  equivalentes  ou  electives  (modules,  vitesses  des  ondes  de 
volume  et  constantes  d  absorption)  ,  Dans  le  cas  d-?  longueurs  d 1  on do  »r ande  s  devant 
la  dimension  des  objets  diffuseurs^  ’  Kuster  et  Toksoz  ont  developpe 

une  raethode  (ne  prenant  pas  en  compte  les  effets  de  diffusion  multiple) 
qui  donne  les  proprietes  des  materiaux  equivalents  en  comparant  les 
amplitudes  monopolaire  ,  dipolaire  et  quadripolaire  exactes  et  effectives. 


Nous  avons  etendu  cette  approche  au  cas  ou  la  longueur  d'onde  devient 
comparable  aux  dimensions  des  ob  je  ts  di  f  fusan  ts  (  cons  ide re  comme  spherique)  ; 
dans  ce  cas,  on  tient  compte  des  resonances  des  particules  qui  causent  des 
resonances  dlargies  dans  les  parane t res  effectifs  des  materiaux.  On  traite 
de  cette  fa^on  les  liquides  contenant  des  bulles,  les  solides  performs,  et 
les  solides  contenant  des  inclusions.  Nos  resultats  sont  verifies  par  les 
mesures  de  Kinra  et  Anand .  De  plus,  des  resultats  anterieurs  concernant  les 
modules  effectifs  de  milieux  composes,  obtenus  dans  la  limite  statique 
(ou  de  basse  frequence),  sont  trouves  ici  comme  cas  particuliers  de  notre 
approche . 

1.  INTRODUCTION. 

We  employ  a  method  first  suggested  by  Ament  [1],  and  later  elaborated  on  and 
applied  to  the  case  of  seismic-wave  propagation  by  Kuster  and  Toks'oz  [2], 
in  order  to  obtain  effective  material  constants  (moduli,  wave  speeds, 
attenuations)  for  microinhomogeneous  media  consisting  of  a  random 
distribution  of  spherical  solid  inclusions,  fluid-filled  cavities  or  gas 
bubbles  in  an  otherwise  homogeneous  solid  or  fluid  host  medium.  Our  method 
is  based  on  the  scattering  of  waves,  and  by  a  comparison  of  monopole,  dipole 
and  quadrupole  scattering  amplitudes  from  the  microinnomogeneities  with  those 
from  a  sphere  of  the  effective  material,  we  were  able  to  obtain  frequency- 
dependent  effective  moduli,  wave  speeds  and  attenuations  which,  in  contrast 
to  previous  work,  contain  the  effects  of  monopole  resonances  of  the  individual 
inclusions,  as  retained  by  us  in  the  long-wavelength  expansion.  Our  theory 
neglects  rescattering  effects  and  is  thus  restricted  to  small  to  moderate 
volume  concentrations  <J>  (<  40%  in  some  experiments).  By  its  dynamic  nature, 
it  furnishes  dispersion  curves  for  wave  propagation  constants  which  reproduce 
well  their  acoustical  branch,  as  obtained  experimentally  for  random  glass 
spheres  in  an  epoxy  matrix  [3],  while  the  optical  branch  (dipole  resonances)  is 
not  included  in  our  results  due  to  our  restriction  to  monopole  resonances. 

We  also  discuss  our  results  for  the  case  of  water  containing  air  bubbles,  as 
well  as  recent  experiments  on  this  topic  [4].  Our  general  dynamic  results 
are  shown  to  reduce  to  previously  obtained  static  limits  for  the  effective 
material  constants. 

2.  GENERAL  FORMALISM. 

In  an  elastic  host  medium  containing  N  solid  spherical  inclusions, 
the  total  eompressional-wave  potential  generated  by  an  incident  p-wave  is 


(2.1) 


(r)  +  ^  ^  , 

U1 

if  one  neglects  multiple  scattering  and  retains  only  single-scattering 
(k)  -*■ 

amplitudes  $sc  (r,  r^)  from  the  N  scatterers  of  radii  a^,  their  centers 
•  •  . 

being  at  random  position^  r^.  Alternately,  one  may  consider  a  large  sphere 
of  radius  R,  with  center  at  r  ,  consisting  of  an  "effective"  medium 
(quantities  denoted  by  a  tilde)  and  imbedded  in  the  same  host  material  ; 
it  leads  to  a  total  field 

fsc  (2.2) 

We  define  the  effective  medium  to  be  such  as  to  produce  the  same  scattered 
far  field  as  the  actual  composite,  so  that  at  r  >>  R, 

N 

T  -*  —  x  ,  (k)  -*  -  . 

fisc  =  2_  fisc  (r,^k).  (2.3) 

k--± 

In  this  latter  relation,  a  comparison  of  monopole  (Aq) ,  dipole  (A^)  and 
quadrupole  (A^)  scattering  amplitudes  (which  can  be  thought  of  as  the 
leading  terras  in  a  long-wavelength  expansion  of  )  is  shown  to  be 
sufficient  to  determine  the  effective  material  constants  [2].  It  is 
necessary,  however,  to  carry  out  a  similar  long-wavelength  expansion  in  the 
known  [5,6]  multipole  amplitudes  A^  (which  would  appear  to  be  called  for 
here,  for  the  sake  of  consistency)  with  a  certain  amount  of  care  :  some 
terms  in  A  contain  the  square  of  the  frequency,  which  is  normally  neglected 
in  the  long-wavelength  expansion,  but  multiplied  by  p  /p^,  the  ratio  of 
host  ( 1 ) -to-inc lus ion  (2)  densities,  which  may  be  large  so  that  the 
corresponding  term  is  no  longer  negligible.  Retaining  such  terms  in  Aq  is 
seen  [7,8]  to  effect  an  inclusion  of  monopole  resonance  terms  in  the 
scattering  amplitude,  which  indeed  appear  prominently  in  such  quantities  as 
the  dynamic  (frequency-dependent)  sound  speeds  and  attenuations  in  bubbly 
liquids  [7],  or  in  the  effective  wave  speeds  and  moduli  of  materials 
containing  fluid-filled  cavities  [C]  or  solid  inclusions  [9].  These  resonance 
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terms  are  retained  here  in  A  ,  but  not  in  A.  or  A0 .  For  the  case  of  solid 

o  i 

inclusions,  this  leads  to  the  explanation  of  the  acoustical  branches  of 
measured  dispersion  curves  [2]  for  the  effective  medium,  but  not  their 
optical  branches  ;  these  would  have  to  be  explained  by  retaining  resonance 
terms  in  A  j . 

Denoting,  for  a  sphere  of  radius  R,  the  normalized  frequency 

x,,  =  k  , ,  R  where  k,,  =  co/c,  is  the  p-wave  number  in  medium  1,  one  has  [8] 
d  i  d  1  d  1  a 


3  ^  2  *  * 

3  ~~  Xof  “  Xos  Xoi 

31  Xdl  U  +  'iXji  )-X0f  -Xos+Xoi 


(2.4a] 


A.-^-U-Sb 


(2.4b] 


a  _ ,  <2.4c) 

1  3/1  ^2.  (k=1+2p-J+  p.1  (9 


where 


X0f  = 


(2.5a) 


is  the  Minnaert  [10]  bubble  resonance  frequency. 


X-s  =  (2<5i/^4  iY 


(2.5b) 


that  of  the  Meyer-Brendel-Tamm  [11]  resonance  for  evacuated  cavities  in 
solids,  and 


xj  = 


(2.5c) 


an  additional  resonance  due  to  shear  in  the  inclusion  ; 
I  "N  .  2 


+  3 


(2.5d) 


is  the  elastic  bulk  modulus  of  medium  1  (X  and  y  being  the  Lame  constants) 
3.  SPECIFIC  RESULTS 

We  may  ncv  compare  (n  =  0,  1,  2)  separately  in  the  fashion  of 
Eq.  (2.3).  This  determines  elastic  moduli  and  wave  speeds  as  listed  in 
Ref.  [9]  .  One  has,  e.g.,  for  the  (normalized)  effective  p-wave  speed  : 


2(U^V") 


[i -a- p,  /p1)i]Vi  v  (xl*  x1)1'*  ( u  W1 ) 
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(3.5c) 


Here,  the  volume  concentration  is 


E- 


TT 


(aDl  )X3LbdAl  , 


(3.6a) 


where  g(a)  is  the  radius  distribution  of  the  inclusions.  From 


9*  -  (l-$  )?i  +  0.7) 

and  Eq.  (3.3),  one  finds  the  effective  s-wave  speed  : 

~  ~  ~  N  4-/Z 

~  (  f^Z  /  <p  2.  )  >  (3.8) 

which  is  frequency  independent. 

Figure  3.1  shows  theoretical  [9]  and  measured  [3]  results  for  the 

frequency  dependence  of  the  effective  p-wave  speed  in  a  medium  of  glass 

spheres  (average  radius  1  mm)  in  epoxy,  at  concentrations  $>  =  5 %  and  15%. 
Figure  3.2a  shows  this  speed  at  f  =  0.5  MHz  plotted  vs.  <J>,  and  Fig.  3.2b 
displays  the  shear  speed  vs .  5  ;  the  agreement  is  good  up  to  $  =  40%. 

4.  PARTICULAR  CASES 

The  static  effective  material  constants  are  obtained  in  the  limit  f  ■+•  0 . 
Our  above  results  then  agree  with  the  previous  static  results  of  Mai  and 
Bose  [12]  ,  Kerner  [13]  , and  Hashin  [14].  For  rigid  inclusions  ( \ ^ »  y9  ■*  ») 

of  dilute  concentration  (<{>  <<1)  we  recover  the  results  of  Moon  and  Mao  [15], 

and  if  the  host  medium  is  additionally  assumed  incompressible  (Poisson  ratio 
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V  2  -j,  one  has 


/Tx/fW  -  *-+f  2\  (4.0 

a  result  obtained  by  Einstein  [16]. 

5.  GAS  BUBBLES  IN  FLUIDS 

Here  one  has  c  .  =  c  „  =  0,  further  U  =  1,  V  =  0.  The  theory  then 
si  s2 

reduces  to  a  previously  treated  case  [7],  and  Fig.  5.1a  presents  numerical 
results  for  water  containing  a  =  1  mm  air  bubbles  at  <J>  -  0.000377  (their 
resonance  frequency  being  3.285  kHz)  compared  to  experimental  results  of 
Silberraan  [17]  :  the  top  portion  shows  absorption,  the  botton  portion  the 
effective  sound  speed.  Figure  5.1b  gives  more  recent  measurements  (crosses) 
at  4>  =  0.0000267  for  a  bubble  distribution  centered  on  a  =  45  ytn,  compared 
to  theory  (solid  curves)  [4],  In  this  latter  work,  it  is  also  shown  how, 
with  one  measurement  of  c  at  low  frequency,  and  two  at  high  frequency,  the 
first  three  moments  of  the  bubble  size  distributions  can  be  experimentally 
determined. 

6.  CONCLUSION 

The  foregoing  theory  permits  the  determination  of  dynamic  effective 
material  constants  for  microinhomogeneous  media  containing  gas  bubbles, 
fluid-fi! ' -  d  cavities  or  solid  inclusions  (treated  as  spherical),  for  low 
to  model  concentrations  of  the  inhomogeneities,  and  taking  their 
resonanc  as  functions  of  frequency)  into  account.  For  higher  concentrations 
sel f-con- . s tent  methods  have  to  be  invoked  such  as  that  of  Chaban  [18]  ;  this 
latter  method,  however,  treats  resonances  in  an  ad-hoc  fashion  only. 

This  work  represents  a  collaboration  of  the  author  with  G.C.  Gaunaurd, 
with  portions  supported  by  NSWC  and  ONR. 
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FIGURE  CAPTIONS 

Fig.  3.1.  Effective  sound  speed  c^  vs.  frequency  (0-2.5  MHz)  for  a  =  1  ram 
glass  spheres  in  epoxy  (concentrations  $  =  5%  and  15%). 

Fig.  3.2.  Effective  wave  speeds  for  a  =  1  ma  glass  spheres  in  epoxy,  vs. 

concentration  $  :  (a)  effective  sound  speed  at  f  =  0.25  MHz, 

(b)  effective  shear  wave  speed. 

Fig.  5.1.  Attenuation  (top)  and  effective  sound  speed  (bottom)  in  bubbly 

water  at  (a)  a  =  1  ram,  $  =  0.000377,  (b)  a  =  45  urn,  $  =  0.0000267. 
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Baum's  "Singularity  Expansion  Method"  (SEM),  formulated  for  radar 
scattering,  is  based  on  the  observation  that  transient  scattered  echoes 
appear  to  be  composed  of  a  sum  of  decaying  sinusoids.  Fourier-trans forming 
these  expressions  into  frequency  space  reveals  a  manifold  of  simple  poles 
in  the  complex  frequency  plane,  identical  with  those  of  the  Resonance 
Scattering  Theory  (RST),  and  commonly  grouped  in  layers.  We  have  carried 
out  a  time  dependent  analysis  of  scattering  from  rigid  and  elastic  spheres 
and  cylinders,  with  the  following  results.  The  use  of  short  pulses  (of 
spatial  extent  small  compared  to  the  scatterer's  dimension)  corresponds  to 
the  echo  being  a  residue  sum  over  a  large  number  of  SEM  poles,  and  the 
residue  sum  over  one  given  pole  layer  produces  a  sequency  of  echo  pulses 
corresponding  to  a  creeping  wave  repeatedly  circumnavigating  the  sphere 
with  the  appropriate  group  velocity.  The  use  of  f ini te- length  sinusoidal 
wave  trains  (long  compared  to  the  scatterer's  extension)  produces  a 
reflected  wave  train,  coherently  superimposed  by  a  sequence  of  overlapping 
creeping  wave  trains,  which  cause  initial  transients  as  well  as  a  final 
transient  tail  following  the  echo.  These  transients  only  appear  if  the 
carrier  frequency  coincides  with  an  eigenfrequency  of  the  target,  and  the 
tail  amplitude  plotted  as  a  function  of  frequency  then  reproduces  the 
spectrum  of  resonances  including  their  widths,  leading  to  a  direct  target 
spectroscopy  as  accomplished  experimentally  by  Ripoche  et  al .  This  tail 
corresponds  to  the  ringing  of  a  given  eigenvibration ,  which  is  selectively 
excited  when  overlapped  by  the  narrow  spectrum  of  the  long  incident  pulse. 


RESUME  :  La  "methode  de  de ve 1 oppement  en  s ingul ari tes  '  ["Singularity  Expansion 

Method"  (SEM)]  de  Baum,  formultle  pour  le  cas  de  la  diffusion  radar,  est 
basee  sur  1 1 observat ion  que  les  echos  de  la  diffusion  transitoire  paraissent 
€tre  composes  d'une  somme  de  sinusoTdes  araorties  .  !.a  trans  format  i  on  de  ces 
sinusoides  dans  l'espace  des  frequences  par  la  transformation  de  Fourier, 
revele  une  multitude  de  poles  sinples  dans  le  plan  cornplexe  des  frequences, 
identiques  a  ceux  de  /.a  "Thdorie  de^  la  diffusion  resonnante"  ["Resonance 
Scattering  Theory"  (RST)  ] ,  et  gener aleraent  groupes  en  couches.  Nous  avons 
execute  une  analyse  de  la  diffusion  transitoire  par  des  spii/res  et  des  cvlindres 
impenetrables  ou  elastiqucs,  avec  les  resultats  suivants.  L1 ut i 1 isat ion 
d' impulsions  courtes  (d'une  extension  spatiale  petite  en  comparaison  avec 
les  dimensions  de  l'objet  diffusant)  correspond  au  fait  que  l'echo  est 
une  somme  sur  les  residus  d'un  grand  nom.bre  de  poles  SEM,  et  la  somme  sur 
les  residus  d'une  couche  de  poles  donnee  produit  une  suite  d' impulsions 
qui  correspond  a  une  ondc  c i rconf erent ie 1 le  encerclant  la  sphere 
plusieurs  tours  avec  la  vitesse  de  groupe  appropriee.  L ' ut i 1 isat ion  de 
trains  d'ondes  sinusoidales  d'une  duree  finie  (d'une  extension  spatiale 
grande  en  comparaison  avec  1'etendue  de  l'objet  diffusant)  produit  un 
train  d'ondes  reflechies,  sur  lequel  est  superposee  d ' une  faqon  coherente 
une  suite  de  trains  d'ondes  circonf erentielles  qui  se  chevauchent  et 

causent  des  transitoires  initiaux  ,  ainsi  ou'  une  trainee  transitoire 
qui  suit  l'echo.  Ces  transitoires  apparaissent  seulement  si  la  frequence 
porteuse  du  train  d'ondes  incident  coincide  avec  une  frequence  propre  de 
la  cible  le  trace  de  la  courbe  donr.ant  1  'amo  1  i  tudo  de  la  trainee  en  fonction  de  la 
frequence  estalors  ident ique a celui  des  resonances,  avec  lour  largeur,  ce  qui 
autorise  une  spectrosconie  directe  de  la  c  ib  le ,  a  ins  i  que  1'ont  montre  Ei  ooeho  e  t  a 
exoerimentalement.  La  trainee  correspond  a  la  reemiss i on  fibre  d 'une  v: b rat  ion 
propre  distincte,  excitee  d'une  faqon  selective  si  sa  frequence  se  trouve 
dans  la  region  etroite  du  spectre  de  1' impulsion  (longue)  incidente. 

1.  INTRODUCTION. 

The  singularity  Expansion  Method  (SEM),  first  formulated  for  the  case 
of  radar  scattering  [1],  is  based  on  the  observation  that  the  echoes  of 
radar  pulses  scattered  from  a  finite  target  appear  as  the  superposition 
of  damped  sinusoids, 

=  (1.1) 
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with  complex  amplitudes  R.  and  exponents  s  .  Taking  the  Laplace  transform, 

Ct  Ct 

this  pulse  shape  then  indicates  the  presence  of  complex-frequency  poles 
in  the  scattering  amplitude. 


(1.2) 


where  s  =  oj/i,  and  u)  =  Im  s  are  the  natural  frequencies  of  the  target 

a  a 

whose  existence  gives  rise  to  resonances  in  the  scattering  amplitude, 
Eq.(1.2).  This  latter  equation  corresponds  to  a  single-frequency  excitation 
of  the  resonances.  If  the  incident  amplitude  is  pulsed  in  time,  it  can  be 
represented  by  the  Fourier  integral 

fhc(  2,-0=  |G(k)eVk("'"t)dl/27r,  (1.3, 

* 

where  k  =  uj/c.  Its  spectrum  G  (k)  weighs  the  factor  in  Eqs.  (1.1),  (1.2). 

A  short  pulse  is  characterized  by  a  wide  spectrum  G(k)  ;  for  G(k)  =  1, 
one  has  f^nc(z,t)  =  <5(z-ct).  Numerous  short-pulse  experiments  for  acoustic 
scattering  from  submerged  elastic  bodies  [2-6],  or  for  elastic-waves 
scattering  from  cavities  [7,8]  have  shown  that  the  echo  then  consists  of  a 
reflected  pulse,  with  a  shape  close  to  that  of  the  incident  pulse  [9], 
followed  by  a  sequence  of  pulses  slightly  widened  as  compared  to  the 
incident  pulse  shape  [10,11],  which  correspond  to  the  multiple  encirclement 
of  the  scatterer  by  surface  (creeping)  waves.  Theoretically,  the  poles 
k  =  is  /c  at  the  scattering  amplitude  fall  into  lavers  in  the  complex 
frequency  plane,  generally  parallel  to  the  real  axis  and  displaced  from  it 
at  successively  increasing  imaginary  spacings.  The  wide  pulse  spectrum 
covers  many  poles  in  the  first  layer  and  in  succeeding  ones,  and  their 
coherent  residue  sums  add  up  to  the  experimentally  observed  time  sequence 
of  pulses  [12],  one  sequence  for  each  layer  sum  which  hence  corresponds 
to  a  given  mode  of  surface  wave.  The  pulses  were  found  to  propagate  with 
the  group  velocity  of  the  surface  waves. 

An  incident  pulse  of  long  duration  has  a  narrow  spectrum  G(k) ,  whose 
weight  in  Eq . (1 .2)  can  radically  limit  the  number  of  poles  contributing 
to  Eq.(l.l).  In  fact,  if  the  pulse  duration  is  chosen  long  enough,  the 
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width  of  G(k)  may  be  made  less  than  the  spacing  between  the  real  parts 
of  the  poles,  =  Re  k^.  If  then  the  peak  of  G(k)  coincides  with  the 
position  k^  of  the  ath  pole,  Eq.(l.l)  becomes 


fsc(D  *  ^  G 


having  been  reduced  to  the  contribution  of  one  single  pole  a,  which 

appears  in  the  form  of  one  single  damped  sinusoid  with  its  decay  factor 

exp(Re  s  )t  =  exp(-  | Im  k  j)ct, Im  k  <0,  representing  the  "ringing"  of 

a  ot  cx  _ | 

the  resonance  a  ;  this  ringing  is  continuous  and  of  duration  (Jin  k  ;c) 

Such  a  use  of  long  pulses  leads  to  the  possibility  of  exciting  target 
resonances  corresponding  to  individual  poles  of  tne  scattering  amplitude, 
and  hence  determining  the  natural  frequencies,  the  imaginary  parts  Im  k^ 
of  the  pole  positions  (from  the  decay  constant  of  the  ringing),  and  the 
residues . 

The  information  thus  obtained  can  conceivably  be  utilized  for  target 
identification  purposes  [13].  The  pioneering  experimen'  of  Maze,  Ripociie  et 
al,  [14-18]  have  clearly  demonstrated  the  ringing  of  target  resonances  as 
excited  by  long  incident  sound  pulses,  thus  obtaining  the  target's 
eigenf requency  spectrum  in  a  direct  fashion,  and  opening  up  the  way  to  what 
may  be  called  "acoustic  spectroscopy"  [13,19], 

2.  SCATTERING  OF  SHORT  PULSES. 

The  acoustic  field  in  the  presence  of  a  sphere  is  [20] 

CO 

P-H; 

y\  =  0 

and  its  poles  are  given  by  those  of  the  S-function 


1=5 


(V) 


US-\ 

1  /n _ 


(2.2) 


(where  x  -  ka,  a  -  sphere  radius)  in  the  complex  frequency  plane.  These 
are  shown  in  Fig.  2.1  [21]  for  the  case  of  a  tungsten  carbide  sphere. 
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Figure  2.1(a),  Cop,  presents  fh^  pules  x  0  corresponding  to  the  resonances 

n-*'  (  i )  f 

of  external  surf-.wc  waves  (close  to  the  series  of  h  (x)  ;  the  pole 

layers  corresponding  to  th e  ,£  th  wave  a- 1,2...  .)  are  connected  by  solid 

lines  (while  the  poles  of  a  given  mode  n  are  connected  by  dashed  curves). 

Part  (b)  of  the  figure,  bottom,  shows  the  poles  due  to  internal  surface  waves 

(given  by  the  zeros  of  F  *  -  z  ^  ,  the  layer  labeled^  =  1 

6  1  n  n 

corresponding  to  the  Rayleigh  wave. 

For  the  case  of  a  rigid  sphere  (F  ->•  0),  one  obtains  by  an  expansion 

( r  )  n 

of  the  denominator  of  around  its  poles  for  the  JL  th  creeping- 

wave  amplitude  in  the  far  field,  after  application  of  Eq.(1.3)  : 


9-rJ4rrr^)T'M)  ’  <2-3) 


-  <x  .r 


where  T  =  (ct-r)/a  is  a  dimensionless  time  variable.  This  is  the  mentioned 
weighted  residue  sum  over  the  SEM  poles.  For  an  incident  5  -  pulse  (G  =  1), 
Fig.  2.2  shows  the  backscattered  pulse  sequence  arising  from  a  numerical 
sum  over  50  decaying  sinusoids  in  Eq.(2.3),  corresponding  to  the  £.  =  1 
external  creeping  wave.  The  sinusoids  are  seen  to  cancel,  except  at  the 
arrival  times  of  the  creeping  wave  which  circumnavigated  the  sphere  over 
1/2,  3/2...  circumferences,  the  arrival  times  being  determined  by  the 
group  velocity  of  the  surface  wave. 

3.  SCATTERING  OF  LONG  PULSES. 

This  technique  was  pioneered  experimentally  by  Maze  and  Ripoche  [14-18], 

Figure  3.1(a)  shows  their  initial  pulse,  a  rectangular  wave  train  long 

enough  to  wrap  around  their  cylindrical  target  50-100  times.  Figure  3.1(b) 

presents  the  reflected  echo  when  the  carrier  frequency  of  the  train 

coincides  with  an  eigen  frequency  ^  =  Rex^  ;  the  figure  shows  both  an 

initial  transient  and  the  pulse's  settling  down  to  a  steady-state  plateau, 

while  after  the  passage  of  the  reflected  train,  a  tail  due  to  the  ringing 

of  the  resonance  appears.  Its  decay  time  determines  Irax  „  ,  and  the  rate 

n  It 

of  disappearance  of  the  tail  as  the  carrier  frequency  is  moved  away  from 
the  resonance  frequency  determines  the  width  of  the  resonance.  In  this  way, 
an  eigenfrequency  spectrum  of  the  cylinder  is  obtained  as  shown  in  Fig.  1.1(b) 
of  Ref.  [13],  by  performing  a  sweep  of  the  carrier  frequency. 


The  physics  of  this  ringing  has  been  studied  recently  both  for 
elastic  spheres  and  cylinders  [22].  Figure  3.2  shows  an  incident  wave 
train  (center),  its  spectrum  overlaid  on,  and  coinciding  with  an  interference 
minimum  of  the  backseat te ring  amplitude  for  a  tungsten  carbide  sphere  (top), 
and  the  backs  cat tered  s ignal  at  resonance  (x^£  =  14.07).  Figure  3.3  shows 
similar  quantities  with  the  carrier  frequency  somewhat  off  resonance  (x  =13,7) 
The  overlay  with  the  minimum  causes  a  deep  constriction  in  the  plateau 
(steady-state  portion)  of  the  reflected  signal  ;  an  initial  transient 
descends  to  the  plateau  and,  after  the  passage  of  the  reflected  pulse,  a 
final  transient  (tail,  =  ringing  of  the  resonance)  decays  in  a  fashion 
similar  to  the  initial  transient,  both  showing  a  pattern  of  steps  first 
pointed  out  by  Veksler  [23], 

The  transients  originate  from  a  coherent  superposition  of  ingredients 
schematically  shown  in  Fig. 3. 4.  Disregarding  the  dotted  rectangles 
representing  internal  reflections  in  backseat te ring ,  one  has  a  specularly 
reflected  wavetrain  (large  rectangle)  and  successive  circumferential 
wavetrains  (small  rectangles).  At  resonance,  all  these  latter  wavetrains 
interfere  constructively  with  each  other  and  form  a  tail  which  consists  of 
steps,  hereby  approximating  the  exponential  ringing,  Eq,(1.4),  of  the 
resonance.  They  also  form  a  coherent  initial  transient,  region  (i)  of 
Fig.  3.4,  which  may  interfere  either  constructively,  destructively  or  in 
an  intermediate  fashion  with  the  specular  echo  depending  on  the  material 
parameters.  In  the  case  of  Fig.  3.2,  the  interference  is  destructive  as 
indicated  by  the  minimum  of  the  form  function  ;  the  interference  of  the 
initial  transient  with  the  specular  echo  thus  leads  to  a  step-wise  descent 
of  the  signal  towards  the  quas i-s tationary  regime  [region  (ii)  of  Fig.  3.4], 
ending  with  a  constricted  plateau  in  Fig.  3.2.  For  the  off-resonance  case 
of  Fig.  3.3,  the  tail  has  decreased  in  amplitude  (indicating  the  resonance 
width),  and  the  interference  of  the  initial  transient  with  the  specular 
pulse  has  changed. 


CONCLUSION. 


A  complete  understanding  of  the  transient  scattering  of  long  or  short 
sound  pulses  from  elastic  targets  is  obtained  by  considering  the  interplay 
of  the  incident  pulse  spectrum  with  the  SEM  pole  pattern  of  target 
eigenfrequencies  in  the  complex  frequency  plane,  and  with  the  backscattering 
amplitude  as  plotted  vs.  frequency.  In  the  time  domain,  each  pole  contributes 
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a  residue  in  the  form  of  a  decaying  sinusoid  -  the  ringing  of  the  correspond!:', 
scattering  resonance.  Short  echo  pulses  consist  of  a  sura  over  many  poles, 
with  their  ringing  all  adding  up  to  nonvanishing  (short)  echo  signals  of 
circumferential  pulses,  with  arrival  times  corresponding  to  the  group 
velocities  of  the  circumferential  waves.  Long  echo  pulses  contain  a  coherent 
sum  of  overlapping  circumferential  wavetrains,  constructively  interfering 
at  resonance  but  destructively  off  resonance,  hence  leading  to  an 
approximately  exponentially  (actually,  step-wise)  decaying  ringing  of  the 
resonance,  or  to  the  absence  of  ringing  off  resonance.  These  trains  interfere 
arbitrarily  with  the  specular-echo  wave  train,  causing  initial  transients  and 
a  steady-state  plateau  in  the  latter.  The  success  of  MI1R  [18]  as  a  tool  of 
acoustic  spectroscopy  [13,19]  rests  in  the  direct  measurement  of  resonance 
amplitudes  (via  the  ringing  effect),  unencumbered  by  any  interference  with 
the  specular-reflection  background  as  in  the  steady-state  regime.  Due  to  the 
prominence  and  sharpness  of  (intrinsic)  elastic-wave  resonances  (i.e.,  their 
high  Q  values),  and  the  ensuing  capability  of  the  associated  surface  waves 
to  encircle  the  target  a  large  number  of  times  before  being  attenuated, 
resonance  spectroscopy  in  underwater  acoustics  may  well  become  a  useful  means 
for  target  recognition  if  applied  properly  and  with  sufficient  insight  and 
care.  ' 

Portions  of  the  work  reported  here  were  supported  by  the  ONR,  NRL, 
and  NSWC . 
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SEM  -  pole  layers  corresponding  to  the  resonances  of  (a) 
external  (top)  and  (b)  internal  (bottom)  surface  waves  on 
a  tungsten  carbide  sphere  in  water. 

Backscattering  pulse  sequence  of  first  {JL=1)  external 
circumferential  wave  on  a  rigid  sphere. 
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4.  Schematic  view  of  specular  (large  solid  rectangle),  penetrating 
and  multiply  internally  reflected  (dotted  rectangles),  and 
circumferential  wave  trains  (small  solid  rectangles),  before 
their  coherent  addition  ;  (i)  initial  transient  region, 

(ii)  quasi  -  steady  state  region,  (iii)  final  transient  region 
(of  ringing) . 

2.  Amplitude  of  incident  wave  train  (center),  its  spectrum  overlaid 

on  the  backscattering  amplitude  of  a  tungsten  carbide  sphere  (top), 
and  backscattered  signal  (bottom),  at  a  resonance  (x^)  =  14.07). 

3.  As  in  Fig.  3.2,  but  off  resonance  (x  =  13.7). 
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SUMMARY  :  Echoes  of  acoustic  or  elastic  waves  scattered  from  a  target  carry 

within  them  the  resonance  features  caused  by  the  excitation  of  the 
eigenvibrations  of  the  target.  By  means  of  a  suitable  background  subtraction 
it  is  possible  to  isolate  the  target's  spectrum  of  resonances.  This  spectrum 
characterizes  the  target  just  as  an  optical  spectrum  characterizes  the 
chemical  element  or  compound  that  emits  it.  Extracting  the  resonance 
information  from  the  echo  allows  the  possibility  of  identifying  the  target 
as  to  its  size,  shape,  and  composition.  This  is  illustrated  by  studying  the 
dependence  of  the  resonance  spectra  of  fluid  targets  upon  changes  of  target 
shape,  including  variations  from  spheres  to  prolate  spheroids  and  finite- 
length  cylinders.  The  resulting  "acoustic  spectroscopy"  (a  concept  introduced 
by  Derem)  generates  the  same  type  of  level  scheme  as  in  optics,  and  it  may  thu 
be  used  for  solving  some  aspects  of  the  "inverse  scattering  problem"  (i.e., 
the  pro’  lem  of  identification  of  the  nature  of  the  target  from  the  returned 
echoes).  A  study  along  these  lines  shows  that  the  eigenfrequencies  of  fluid- 
filled  cavities  in  a  solid  medium  (obtained  by  us  in  the  complex  frequency 
plane)  are  tied  to  resonance  features  in  the  scattering  amplitude  which  can 
be  analyzed  to  provide  the  material  parameters  (density  and  sound  speed)  of 
the  fluid  filler  :  the  resonance  spacing  giving  the  sound  speed,  and  the 
resonance  widths  the  fluid  dens ity-hence  leading  to  a  solution  of  the  inverse 
scattering  problem  in  this  case. 

RESUME  :  Les  echos  d'ondes  acoustiques  ou  elastiques  diffuses  par  une  cible 

portent  en  eux  des  effets  de  resonances  causes  par  l'excitation  des 
vibrations  propres  de  la  cible.  Par  le  moyen  d ' une  soustraction  appropnee 
du  fond,  il  est  possible  d'isoler  le  spectre  des  resonances  de  la  cible.  Un 
tel  spectre  caracterise  la  cible  aut  int  qu'un  spectre  optique  caracte  rise- 
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l'dlement  ou  le  compose  chimique  emetteur  de  ce  spectre.  L'extraction  des 
echos  de  1 ' in f orma t i on  sur  les  resonances  offre  la  possibility  d' identifier 
la  cible,  par  sa  grandeur,  sa  forme,  et  sa  composition.  Nous  illustrons  cela 
en  etudiant  la  dependence  des  spectres  de  resonances  de  cibles  fluides  en 
changeant  la  forme  de  la  cible,  ceci  inclus  les  variations  entre  une  spi.ere 
et  un  spheroi'de  oblong ,  ou  un  cylindre  de  longueur  finie.  Le  concept  a  suivre 
de  "spectroscopie  acoustjque"  (concept  introduit  par  Derem)  a  comine  base  la 
meme  sorte  de  schema  de  niveaux  que  1 1  on  trouve  en  optique,  et  il  peut 
alors  etre  utilise  pour  resoudre  en  quelque  sorte  le  "probleme  inverse" 

( c ' es t- a- dire  1 1  identification  de  la  nature  de  la  cible  en  etudiant  les  echos). 
Une  telle  etude  montre  que  les  frequences  propres  des  cavites  dans  un  milieu 
solide  remplies  d ' un  fluide  (obtenues  ici  dans  le  plan  complexe)  sont  liees 
aux  effets  de  resonances  dans  1 ' amplitude  du  signal  diffuse,  qui  peuvent  etre 
analysees  afin  de  fournir  les  parametres  physiques  du  fluide  :  l'espacement 
des  resonances  donnant  la  vitesse  du  son,  la  largeur  des  resonances  donnant 
la  densite  du  fluide,  menant  ainsi  a  la  solution  du  probleme  inverse  dans  ce 
cas  . 

1.  INTRODUCTION. 

In  1979,  A.  Derem  wrote  [1]  :  "  L'analyse  de  la  diffusion  des  ondes 
acoustiques  par  un  cylindre  elastique  immerge,  fait  surgir  une  veritable 
spectroscopie  acoustique".  This  was  said  following  an  analysis  of  effects  of 
the  resonances  of  elastic  objects,  after  their  excitation  by  incident  sound, 
as  they  appear  in  the  scattered  echoes,  and  it  seems  to  be  the  first 
formulation  of  this  concept  in  conjunction  with  elastic  resonance  frequency 
schemes.  (The  previous  terra  "ultrasonic  spectroscopy",  and  the  related  early 
experiments  of  Gericke  [2,3],  concern  the  more  general  notion  of  the 
ultrasonic  spectrum  reflected  from  scatterers  in  solids,  and  the  information 
on  these  scatterers  contained  in  the  spectrum,  not  referring  specifically 
to  any  resonances).  The  new  concept  was  brilliantly  justified  experimentally 
shortly  afterwards  by  51a ze  ,  Ripoche  and  others  [4,5]  in  their  pioneering 
experiments  which  exhibited  in  a  direct  fashion  the  mechanical  eigenvibrations 
of  a  submerged  elastic  cylinder,  together  with  the  "level  scheme"  of  the 
corresponding  resonances  excited  by  incident  sound,  straightforwardly  extracted 
from  the  reflected  echo  [6],  Figure  1.1  shows  the  resonances  in  the  reflected 
sound  of  an  aluminum  cylinder  in  water,  as  contained  in  the  total  echo  (top), 
and  as  extracted  e xpe r imen ta 1 1 v  (bottom)  [7].  The  bottom  curve  represents 
the  "acoustic  spectrum"  of  the  cylinder,  which  appears  strikingly  similar 
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to  the  optical  spectrum,  or  level  scheme,  of  atoms  or  chemical  compounds. 

These  facts  then  raise  the  question  whether,  as  in  the  optical  case, 
such  a  measured  spectrum  can  provide  information  on  the  source  of  the 
spectrum  (here,  the  scattering  object).  This  led  to  a  study  of  the 
sensitivity  of  the  acoustic  level  scheme  for  internal  [3,9]  or  external 
resonances  [10,11]  to  changes  in  the  shape  of  the  target.  As  to  the  informat 
which  can  be  extracted  from  the  spectrum  regarding  tne  cons'stency  of  the 
target,  we  have  carried  out  illustrative  calculations  on  fluid-filled 
spherical  [12,13]  or  cylindrical  [14]  cavities  in  a  solid,  or  on  fluid 
cylinders  in  a  fluid  [15],  which  showed  that  in  general,  for  a  cavity  of 
known  shape,  the  sound  velocity  in  the  filler  fluid  can  be  obtained  from  the 
frequency  values  of  the  resonances  (  together  with  the  size  of  the  cavity), 
and  the  density  of  the  filler  fluid  from  tne  width  of  the  resonances. 
Finally,  the  determination  of  the  physical  properties  of  ocean  bottom  layers 
has  also  been  shown  to  be  obtainable  i 1 6  —  1 S  ]  from  the  effects  of  lover 
resonances  contained  in  the  acoustic  echo. 

The  proposed  solution  of  inverse  scattering  problems  by  utilizing 
target  resonances  represents,  of  course,  only  one  particular  approach 
towards  the  solution  of  inverse  problems  [ 1 9  J  .  Other  notable  techniques  in 
ultrasonics  are  the  mentioned  use  of  ultrasonic  spectroscopy  [20],  or  3 
method  of  analysis  to  determine  the  geometry  of  material  defects  based  on 
the  Kirchhoff  approximation  [21,22] . 

2.  ACOUSTIC  RESONANCE  LEVEL  SCHEME. 

A  calculation  of  the  eigenfrequencies  for  prolate  fluid  spheroids  and 
finite-length  fluid  cylinders  in  vacuum,  which  will  approximate  the  internal 
resonance  frequencies  in  another,  much  less  dense  fluid  (except  for  the 
imaginary  parts  of  the  frequencies,  which  vanish  for  ambient  vacuum)  was 
carried  out  [8,9]  in  order  to  gauge  the  sensitivity  of  the  level  scheme 
under  charges  in  shape.  For  the  spheroids,  the  calculation  was  carried  out 
by  subjecting  spheroidal  wave  function  to  the  appropriate  (Dirichlet) 
boundary  condition  on  the  inner  surface  ;  for  f ini te- length  cylinders,  exact 
solutions  are  available.  Figure  2.1  shows  the  e i gen f re quencv  "level  scheme" 
for  spheres  and  infinite  cylinders  as  limiting  cases,  as  well  as  finite- 
length  cylinders  and  spheroids  (axes  ratios  b/a  =  1.11... and  4.0).  The 
striking  "splitting"  (lifting  of  degeneracy)  of  the  limiting-case  levels  fo- 
the  smaller  symmetry  of  finite  non-spher ical  bodies,  the  azimuthal  quantum 
number  no  layer  being  degenerate  as  for  the  sphere  case,  was  explained 
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later  [11,23]  in  terms  of  the  formation  of  helical  surface  waves  with 
discrete  pitch  angles. 

For  the  case  of  penetrable  objects  in  a  fluid  or  solid  medium  (including 
soft  or  rigid  objects  where  only  external  resonances  appear),  one  may  extend 
the  level  scheme  of  Fig.  2.1.  to  a  plot  of  resonance  frequencies  in  the 
complex  frequency  plane  ;  see,  e.g.,  Figs.  4.1  or  5.1  of  Ref.  [23]. 

3.  EXAMPLES  OF  INVERSE  SCATTERING. 


It  has  been  shown  analytically  [12,13]  for  the  example  of  a  fluid-filled 
spherical  cavity  in  a  solid  that  a  determination  of  the  resonance  frequencies 
can  determine  the  sound  velocity  in  the  fluid,  and  of  the  resonance  width 
the  density  of  the  fluid.  Numerical  calculations  for  a  fluid-filler  cylinder 
in  steel  ^14j  bear  this  out.  Figure  3.1  shows  the  dependence  of  the  reduced 
resonance  frequency  spacing  Ax^£  (x  =  ka,  k  =  wave  number  in  the  fluid, 
a  =  cavity  radius)  on  the  ratio  c^/ c2  °f  internal  sound  speed  c^  to  external 
p-wave  speed  c 2>  for  different  values  of  intemal-to-external  density  ratio, 
n  being  the  mode  number  and  X  the  order  of  the  resonance  frequency  within 
the  n  th  mode.  The  near-linear  dependence  on  c^/e2>  an<*  intens^t^v^ty 
to  Pj/p2  indicates  the  feasibility  of  determining  Cj/c2  from  the  resonance 
spacing. 

In  Fig.  3.2,  we  plot  the  resonance  width  Tnn  of  the  n  =  0  resonances 
2  2 

vs,  Pj  Cj,  /p2  c 2  ,  for  different  values  of  This  illustrates  the 

feasibility  of  determining  p^/p2  from  the  resonance  widths.  v 

4.  CONCLUSION 


The  foregoing  represents  an  elaboration  on  the  concept  of  Acoustic 
Spectroscopy,  first  pronounced  by  Derem  [1],  as  a  means  of  determining  the 
properties  of  acoustic  targets  from  the  information  on  the  target  resonances 
which  is  contained  in  the  observed  echo  returns.  The  feasibility  of  this 
concept,  although  yet  to  be  generally  demonstrated  for  cases  of  practical 
interest  [24],  appears  evident  for  targets  of  large  impedance  contrast  with 
the  environment  (so  that  pronounced  resonances  result).  It  is  also  obvious, 
though,  that  an  application  of  this  approach  has  to  proceed,  at  this  time, 
with  complete  understanding  of  the  physics  of  the  scattering  process.  If  this 
caution  is  disregarded,  failure  may  easily  occur.  After  sufficient  development 
of  the  approach,  one  may  ultimately  expect,  however,  that  automated  systems 
for  experimental  acoustic  resonance  spectroscopy  will  have  been  devised. 

Portions  of  this  work  were  supported  by  the  Office  of  Naval  Research, 
the  Naval  Surface  Weapons  Center,  and  the  Army  Research  office. 
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FIGURE  CAPTIONS. 


Fig.  1.1.  Acoustic  backscattering  amplitude  (total  :  top  ;  resonant  : 

bottom)  vs.  frequency  for  a  solid  aluminum  cylinder  in  water. 

Fig.  2.1.  Eigenfrequency  "level  scheme"  of  fluid  spheres,  cylinders 
and  spheroids  in  vacuum. 

Fig.  3.1.  Spacing  of  the  eigenf requencies  of  a  cylindrical  fluid-filled 

cavity  in  steel,  vs.  p-wave  speed  ratio,  for  various  density  ratios 

Fig.  3.2.  Resonance  widths  of  a  cylindrical  fluid-filled  cavity  in  steel, 
vs.  pc  ratio,  for  various  p-wave  speed  ratios. 


